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Abstract: We present a detailed analysis of the classical Dicke-Jaynes-Cummings-Gaudin inte- 
grable model, which describes a system of n spins coupled to a single harmonic oscillator. We focus 
♦pH, on the singularities of the vector-valued moment map whose components are the n + 1 mutually 

,JZj ' commuting conserved Hamiltonians. The level sets of the moment map corresponding to singular 

values may be viewed as degenerate and often singular Arnold-hiouville torii. A particularly in- 
teresting example of singularity corresponds to unstable equilibrium points where the rank of the 
moment map is zero, or singular lines where the rank is one. The corresponding level sets can be 
described as a reunion of smooth strata of various dimensions. Using the Lax representation, the 
associated spectral curve and the separated variables, we show how to construct explicitely these 
level sets. A main difficulty in this task is to select, among possible complex solutions, the physically 
admissible family for which all the spin components are real. We obtain explicit solutions to this 
problem in the rank zero and one cases. Remarkably this corresponds exactly to solutions obtained 

C^) . previously by Yuzbashyan and whose geometrical meaning is therefore revealed. These solutions 

can be described as multi-mode solitons which can live on strata of arbitrary large dimension. In 
these solitons, the energy initially stored in some excited spins (or atoms) is transferred at finite 
times to the oscillator mode (photon) and eventually comes back into the spin subsystem. But 
their multi-mode character is reflected by a large diversity in their shape, which is controlled by 
the choice of the initial condition on the stratum. 



1 Introduction. 

The Dicke model has been used for more than fifty years in atomic physics to describe the inter- 
action of an ensemble of two- level atoms with the quantized electromagnetic field pp. Recently, 
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many experimental and theoretical works have considered systems in an optical cavity where mat- 
ter is coupled to a single eigenmode of the field [21 [3] . These systems offer some unprecedented 
opportunities to realize quantum bits and to process quantum information. The integrable version 
of the monomode Dicke model considered in this paper [H [5] is obtained by making the so called 
rotating wave approximation in which the non resonant terms (e.g. for which a photon is created 
while an atom is excited) are discarded. In the case of a single atom, this model has been solved 
by Jaynes-Cummings [6], and has been used to study the effect of field quantization on Rabi os- 
cillations [7]. In the many atom case, besides the development of cavity QED [2] and circuit QED 
[3] already mentioned, a surge of interest in the Dicke-Jaynes-Cummings-Gaudin model has been 
motivated in the context of cold atoms systems, where a sudden change in the interactions between 
atoms is achieved by sweeping the external magnetic field through a Feshbach resonance [SJ [9] . The 
subsequent dynamics of the atomic system is well described by the DJCG model prepared in the 
quantum counterpart of a classically unstable equilibrium state. 




Figure 1: The image of the moment map for the two spins model. The green dots refer to stable 
(elliptic) equilibrium points. The red dots are unstable (focus-focus) equilibrium points. The 
preimage of a point on a green line is a Liouville torus degenerated to a circle on which the rank 
of the moment map is one. The preimage of a point on a red line is a pinched torus times a circle. 
The parameters (61,62) are chosen so that we have two elliptic and two focus-focus critical points. 



This motivates the study of unstable equilibrium points in classical and quantum integrable 
models. From the classical viewpoint, an integrable system with n degrees of freedom is character- 
ized by a collection of n mutually commuting functions Hi,H2,.-.,H n over the 2n-dimensional phase- 
space. Usually, the physical Hamiltonian H is expressed as a function of these conserved quantities. 
Besides the classical trajectories generated by H, it is of great interest to consider the so-called 
moment map, which associates to any point p in phase space the n-vector (Hi(p), H2(p), ..., H n (p)). 
The image of the moment map is of great importance. Fig.pQ provides an example of such an 
image of the moment map for the two spins Dicke-Jaynes-Cummings-Gaudin model. To each point 
(Hi, H2, ■ ■ ■ , H ri ) in this image corresponds a level set in phase space, that is to say the set of points 
p such that (Hi(p) = Hi,H2(p) = H2, ..., H n (p) = H n ). For many systems of physical relevance 
these level sets are compact (as it is the case here). In the case of a regular value of the moment 
map, the Arnold-Liouville theorem states that the corresponding level set is an n-dimensional torus. 
Non regular values of the moment map, (for which the rank of the moment map may be < n), 
are in some sense more interesting, because they correspond to special torii whose dimension may 
be lower than n or which may contain singular submanifolds. The image of the moment map and 
its critical values encode all the information on the fibration of phase space in Liouville torii. It 
is important to realize that the level set of a non-regular value of the moment map is in general 
a stratified manifold on which the rank of the moment map may jump. The dimension of each 
stratum is equal to the rank of the differential of the moment map on that stratum. 

An extreme case is given by the level set containing an equilibrium point. By definition, the 
gradient of the moment map vanishes at such point (the rank is zero), which is then an equilibrium 
point for any Hamiltonian expressed in terms of Hi, H2,..-,H n . For a purely elliptic equilibrium 
point, the level set is reduced to this point. But for an unstable point, the level set contains the 
equilibrium point, of course, but also a whole manifold whose description is far from trivial. In 
the vicinity of the unstable critical point, the stratification can be qualitatively understood by the 
quadratic normal form [10], deduced from the Taylor expansion of the moment map to second 
order in the deviations from the critical point. A detailed derivation of these normal forms in the 
Dicke-Jaynes-Cummings-Gaudin model has been presented in our previous contribution [11]. But 
the information contained in the normal form is purely local, and it is not sufficient to describe 
globally the level set containing an unstable critical point. Global characteristics can be access by 
viewing these level sets as collections of trajectories, solutions of the equations of motion. In this 
paper, we shall present explicit solutions of the physical Hamiltonian equations of motion for an 
arbitrary initial condition chosen on such a critical level set. They can be described as solitonic 
pulses in which the energy initially stored in some excited spins (or atoms) is transferred at finite 
times to the oscillator mode (photon) and eventually comes back into the spin subsystem. 

We achieve this by using the well known algebro-geometric solution |13l [T5l Q3] of the classical 
Dicke-Jaynes-Cummings-Gaudin model. A key feature of this approach is that the knowledge of 
the conserved quantities Hi,H2,-..,H n+ i, (equivalently a point in the image of the moment map) is 
encoded in a Riemann surface, called the spectral curve. It is of genus g = n, where n is the number 
of spins. The state of the system can be represented by a collection of g points on the spectral 



curve. Remarkably, the image of this collection of g points by the Abel map follows a straight line 
with a constant velocity on the Jacobian torus associated to the spectral curve. This construction 
linearizes the Hamiltonian flows of the model. 

It is well known that critical values of the moment map correspond to degenerate spectral 
curves |16( [13], and this seems to be the most efficient way to locate the critical values of the 
moment map. This strategy has been used extensively in our previous discussion of the moment 
map in the classical Dicke-Jaynes-Cummings-Gaudin model [11] . We believe that the separated 
variables, which we will use in their Sklyanin presentation [15], are a powerful tool to uncover the 
stratification of a critical level set. In this paper, we show that each stratum of dimension n+l — riQ 
on such level set is obtained by freezing uq separated variables on some zeros of the polynomial 
which defines the spectral curve. 

We will aslo pay a special attention to the problem of finding the physical slice on each stratum. 
This is not a straightforward question, because the algebro-geometric method is naturally suited 
to solve equations of motion in the complexified phase-space. But only solutions for which all spin 
components and the two oscillator coordinates are real are admissible physically. The determination 
of the physical slice for a generic regular value of the moment map remains an open problem. But 
on the critical level sets of the moment map, corresponding to a spectral curve of genus zero, the 
problem greatly simplifies. This corresponds to level sets of the moment map containing critical 
points of rank zero or one for which we are able to derive an explicit solution for an arbitrary 
number of spins. Remarkably these solutions were already constructed by Yuzbashyan [12], and 
they find here their place in a broader picture aiming at organizing all particular solutions of the 
Dicke-Jaynes-Cummings-Gaudin model. 

This paper is organized as follows. Section [2] presents the classical Dicke-Jaynes-Cummings- 
Gaudin model of n spins coupled to a single harmonic oscilator, together with the main features of 
the algebro-geometric solution. In particular, we emphasize that the separated variables capture 
in a very natural way the stratification of a critical level set, each stratum being characterized by 
a subset of frozen variables on some double roots of the spectral polynomial. These notions are 
then illustrated in section [3] on the example of the system with only one spin. In this case, the 
reality conditions can be worked out explicitely, both on the critical level sets associated to un- 
stable equilibrium points and for generic values of the moment map. Section [5] gives the complete 
determination of the real slice for the level set of an equilibrium point, and for arbitrary values of n. 
This corresponds to the case where all the roots of the spectral polynomial are doubly degenerate. 
We recover here the normal solitons first constructed by Yuzbashyan [12] . Section [5] addresses 
the same question, in the slightly more common situation where the spectral polynomial has n 
doubly degenerate roots and two simple roots. In this case, we are describing a level set whose 
smallest stratum has dimension 1. The corresponding values of the moment map lie on curves in its 
n-dimensional target space. We benefit in this case from the fact that the spectral curve remains 
rational. We get therefore a general formula for the anomalous solitons in Yuzbashyan's terminol- 
ogy. Section [6] illustrates this general theory for the system with two spins, and section [7] for the 
system with three spins. We present explicit examples of solitonic pulses corresponding to solutions 



living in strata of dimensions 2, 3, and 4. The presence of several directions of unstability away 
from the corresponding equilibrium points is reflected by the non-monotonous time dependence 
seen in these pulses, which contrasts with the simple shape observed for the model with a single 
spin. Our conclusions are stated in section [51 Finally, an Appendix provides a brief summary of 
the construction [TT] of quadratic normal forms of the moment map in the vicinity of equilibrium 
points. 

2 The classical Dicke-Jaynes-Cummings-Gaudin model. 

2.1 Integrability 

This model, describes a collection of n spins coupled to a single harmonic oscillator. It derives from 
the Hamiltonian: 

n n 

H = Y,( 2e J + w ) s i + <*M> + 52 { h J + bs t) (*) 

3=1 3=1 

The Sj are spin variables, and b, b is a harmonic oscillator. The Poisson brackets read: 

{ s -,s b 3 } = -e abc s c 3 , {b,b} = i (2) 

The Sj brackets are degenerate. We fix the value of the Casimir functions 

s j ' s j = s 

Phase space has dimension 2(n + 1). In the Hamiltonian we have used s- = sj ± is 2 - which have 
Poisson brackets {s|, s • } = ±is • , {s^", s^"} = 2is 2 -. The equations of motion read: 

(3) 

(4) 
sf i(2ej+oj)sj -libs] (5) 

■ s - = -i(2ej + u)sj + 2ibs z i (6) 



b = 


= — iub - 


n 
3=1 


S 3 = 


-- i(bsj- 


~ bsf) 



Integrability is revealed after introducing the Lax matrices: 



L(A) = 2Xa z + 2(ba + + ba-) + ^2^-^- (7) 

3=1 J 

M{\) = -i\a z -i^-a z - i{ba + + ba~) (8) 



where a a are the Pauli matrices, a = ^(ff 1 ± io y ) 
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It is not difficult to check that the equations of motion are equivalent to the Lax equation: 



Letting 



we have: 



L(A) = [M(A),L(A)] 



M(A) B(X)\ 
L{X) -{C(X) -A(X)J 



A(X) = 2A + £-^- 

j = l J 



B(X) = 2b + Yl 



(9) 



A — e,: 

sT 



C(X) = 21 + J2tz 



i=i 



A — e,- 



(10) 

(11) 
(12) 



These generating functions have the simple Poisson brackets: 

" Pl2 



{L 1 (X),L 2 ( f x)} = -i 



X — fjL 



,L 1 (X)+L 2 (fi) 



where 
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It follows immediately that Tr (L 2 (A)) = 2A 2 (A) + 2B(A)C(A) Poisson commute for different values 
of the spectral parameter: 

{Tr(L 2 (A 1 )),Tr(L 2 (A 2 ))} = 



Hence 



A(A)^Tr(L 2 (A)) = A 2 (A) + B(A)C(A) 



generates Poisson commuting quantities. One has 

Q2n+2(A) 



A(A) 



71 R; A S 2 



n 7 (A- ei ; 



4A 2 + 4H n+1 +2J2 — -j- + J2 



i A -^- U {x - e > )2 



(13) 



where Q2n+i(A) is a polynomial of degree 2n + 2. The (n + 1) commuting Hamiltonians Hj, 
j = 1, • • • ,n + 1 read: 



)i 



iT n+ i = &6 + J>! (14) 



i=i 



and 



fl i = 26 i ^ + (6*+ + 6 a T) + ^^-^L, i = l,...,n (15) 

The physically interesting Hamiltonian eq.([T]) is: 



*# Cj " €k 



H = uH n+1 + Y,H j (16) 

3=1 

On the physical phase space the complex conjugation acts as s z - = s z and s^ = s~, and of 
course b is the complex conjugate of b. Hence for A real, A(X) is real and -B(A) = C(A). It follows 
that on the physical phase space, Q2n+2(A) is positive when A is real. 

2.2 Separated variables. 

The Lax form eq. ([9]) of the equation of motion implies that the so-called spectral curve T, defined 
by det(L(A) — j£) = is a constant of motion. Specifically: 



r : ^ - A 2 (\) - B(X)C(\) = 0, i.e. M 2 = ^f, \i ( 17 ) 



e i, 



Defining y = /j, Y\j(^—£j), the equation of the curve becomes y 2 = Q2n+2(A) which is an hyperelliptic 
curve. Since the polynomial Q2n+2(A) has degree 2n + 2, the genus of the curve in n. Because 
the model is integrable, it is possible to construct action-angle coordinates (at least locally), but 
their connection to the initial physical dynamical variables is rather complicated. In this work, 
we prefered to work with the so-called separated variables which have the double advantage that 
their equations of motion are much simpler than the original ones and that they are not too far 
from the physical spin and oscillator coordinates. The separated variables are g = n points on the 
curve whose coordinates {Xkil^k) can be taken as coordinates on phase space. They are defined as 
follows. Let us write 

C(X) =2b + ±-^- = 26 nLi(A-A fc ) 
~i X ~ e J n j= i(A-ej) 

the separated coordinates are the collection {\k,^k = A(\k)). They have canonical Poisson brackets 

{Afc,/ifc'} = -iSk,k> 



Notice however that if b = theses variables are not well defined. 

There are only 2n such coordinates which turn out to be invariant under the global U(l) rotation 
generated by H n+ \: 



b^e ie b : b^e~ ie b, sj^e i9 sj, af -»• e~ i6 a 



.1 



So they describe the reduced model obtained by fixing the value of H n+ \ and taking into consid- 
eration only the dynamical variables invariant under this U{1) action. The initial dynamical model 
can be recovered by adding the phase of the oscillator coordinates b, b to the separated variables. 
The equations of motion with respect to H{ in these new variables are (no summation over k): 



OtM - Wk{o )ik, &kj - t — , (o )jp - Yi — r c TTTf — T\ P> 



(19) 



Using the identities: 



c>: 



yZn — r* ; = S P,n-i, o< P <n-i 



the flow associated to the physical Hamiltonian eq. ([T6l) reads: 



9Ak = ^fc 1 - r 77 TT, Hk = A(\ k ) (20) 

Ih^kAk - M) 



The equations of motion eq. (|20p must be complemented by the equation of motion for b, b which 
allows to recover the motion of the full model from the motion of the separated variables in the 
reduced model. We have: 

d t b = iu;b + i Y^ C(X)=iojb-i ^ C(A) = ib(u - 2(Ei - <n(e))) (21) 

Res €j Res oo 

where Si = ^ Aj and <7i(e) = ^ e^. Of course, we also have the complex conjugated equation of 
motion for b. 

We show now how to reconstruct spin coordinates from the separated variables {Xki^k) an d 
the oscillator coordinates (b,b). For C(A) we have eq. (|18p . It is a rational fraction of A which has 
simple poles at A = ej whose residue is s^ . For A(X), we write: 

A(X) - P " +l(A) 



flLi&-ej) 



The polynomial P n+ \(X) — 2A n?=i(A ~~ e j) ^ s °f degree n — 1, and we know its value at the n points 
Xj because: 

n 

P n+1 (X J ) = fi J Yl(\ j -e k ) (22) 

fc=i 
Therefore, we can write: 



n n Y\ ( X — X ) 

P n+1 (X) = 2X JJ(A - ej) + £> " 2A.) J]^ - ^ ff^Ai - A t ) 



(23) 



j=l i fc=l 



Once A(A) and C(A) are known, we can find the spin themselves by taking the residues at A = e,\ 
We get: 



, + _oi iir = ifo-AO ^_ nr=ife-AQ v 



(2A^) n ^ 



n^i(Ai - A/) 



(24) 



The modulus of b is invariant under the global U(l) action generated by H n+ \ and is obtained from 
(recall that in the reduced model H n+ \ is a parameter) 

lb = H n+l - ^ s j 
j 

The phase of b and b is not determined in the reduced model. The same phase appears in the 
formula for s^. Finally, the sj components are obtained from the constraints: 



[st ) 2 + 4s~ = s 2 



They are determined up to the phase of 1/6. 

Note that if the Aj's are arbitrary complex numbers, sj is not the complex conjugate of s^~, 
as it is the case for the physical phase-space. This shows that separated variables are natural 
coordinates on the complexified phase-space. In these variables, the problem of identifying the sets 
{Aj, Hi} corresponding to physical configurations (i.e. those for which all spin components are real) 
is rather non-trivial, and most of the present work is dedicated to it. 

In many situations, and in particular in this article, we are interested in the system with 
prescribed real values of the conserved quantities, Hi,...,H n+ i, i.e. we take as coordinates the 
(\j,Hj) instead of the (Xj,fj,j). This amounts to fixing the Liouville torus we work with, or 
equivalently the spectral polynomial Q2n+2(A). In this setting, the fij are determined by the 
equation 



~[(Xj - e k )nj = ±JQ 2 n+2(Xj) 



fe=l 



The choice of the signs in this formula is very important and will be a recurrent theme in our 
subsequent discussions. They affect the formulae for s z , or equivalently the polynomial P n+ i(A). 
The values of sj's can then be obtained from -B(A) which is determined by: 

, Q 2w+ 2(A)-P n 2 +1 (A) 1 Q 2w+2 (A)-P n 2 +1 (A) 

{) C(\)Y\ k (\-e k )* 26n t (A-A,)n fc (A-6 fc ) {) 

The polynomial in the numerator is of degree 2re, and moreover it is divisible by Yl^X — Aj) because 
P% +1 (k) = <?2n+2(Ai), so we can write 

B(X) = ^f~ Xi \ (26) 

llfel A - e k) 

The spin components «7 are obtained by taking the residue at A = e,- : 

s - = 2b Mi ~ AQ 

3 rifc^'fe'-efc) 

In this consruction, the variables Aj and the corresponding Jii are complicated functions of the 
(Ai, Hi) and of the choice of signs used in the determination of the variables fi{. The set of physical 
configurations (often refered to here as the real slice) is obtained by writing that the set (Aj,/Zj) is 
the complex conjugate of the set (Aj, Hi). This leads to a set of complicated relations whose solution 
is not known in general. It is the purpose of this work to show how to implement them in particular 
cases. 

2.3 Moment map and degenerate spectral curves 

In this section, we discuss the case of a spectral polynomial having double roots, so that Q2n+2(A) 
can be written as: 

n—m 

Q 2n+2 (A) = 4 p 2m+2 (A) JJ (A - Ei)\ P2m+ 2(X) = A 2m+2 + b 2m+1 X 2m+1 + ■■■ (27) 

i=i 

The Ei are either real or come in complex conjugated pairs. Moreover p 2m + 2 (A) is positive for real 
A. As discussed in some previous works [5lUHfT6]. this corresponds to critical values of the moment 
map. The preimage of the moment map for such value is a stratified manifold, whose various strata 
have dimensions ranging from m + 1 to n + 1. Separated variables provide a very natural access 
to this stratification, because each stratum of dimension n — no + 1 can be obtained by freezing 
no <n — m separated variables A^ on no double zeroes of Q 2n + 2 (A). On such stratum, the rank of 
the moment map also drops to n — n-o + 1- As explained before [5], the dynamics of the system for 
initial conditions lying in this stratum is very similar to the one of an effective model with n — no 
spins, but we shall not use this physically appealing result here. 

10 



To produce a polynomial Q2n+2(A) of the form eq. ([2T|) is not completely straightforward because 
it must also be compatible with eq. (|13p . Let us explain how it works. 

The fact that the rational fraction Q2T1+2W /Ylj{^ — £j) 2 has double poles at A — ej whose 
weight is s 2 imposes n conditions: 

n(A - ^)i CJ = * a3 " ife--**) 

. <l\/P2m+2{ej) k¥=j 

where ctj = ±1. The values of the polynomial JL(A — Ei) at the n points ej are therefore known, 
and we can write by Lagrange interpolation formula: 

But the degree of the left hand side is n — m while the degree on the right hand side is superficially 
n — 1, so we have the consistency conditions: 

V — . 3 J = 0, k = 0,---,m-2 (29) 

j 2 v / p 2m+ 2(ej) 

and: 

v S -^£L =1 (30) 

i 2y/p 2 m+2(e j ) 



These equations are obtained by writing that: 

2 ™ U7=i( z - e i) 



Sk,m-x, k = 0, ■■■,m-l 



where the integrals are taken along a contour in the complex plane which encircles all the e,'s. These 
are exactly m conditions on the 2m + 3 coefficients of J>2m+2(A), of which the leading coefficient 
is known. One more condition is obtained by writing that in the right-hand side of eq. (|13p the 
coefficient of A vanishes. This gives: 

b2m+i + 2(ai(e)-a 1 (E)) = 

But we also have: 

dz z m n£T(* - E i) _ ,, r ^ _ V- sa ^ m 



/ 2 ™ Y\t=i\ z -£i) , 



^j 



2y / P2m+2{ej) 
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where the integral is again on a contour encircling all the ey's. Assembling the previous two 
equations gives the consistency condition: 



sa^e 

?2m+l 

3 



*— \ SO! ■)■€„■ 

+ E / J , , =0 ( 31 ) 



Altogether, we are left with m + 1 free coefficients in P2m+2(A) which play the role of conserved 
Hamiltonians of an effective system with m spins [5]. All the other symmetric functions ak(E), 
k > 1 are then determined by eq. ([28l) . 



If m = —1, that is to say j>2m+2(A) = 1> eq. (|28p must be slightly modified as 

[](A - Ek) = (A + 0) H(X - Ci ) + E ^ II( A " e *) ( 32 ) 

* 7 3 i¥=j 

This means 

Comparing with eq. (|13|) at A = oo we see that (3 = 0. 

Let us discuss now the freezing of the A& at the double roots of Q2n+2(A). The equations of 
motion eg. ([20]) become 

„. ,.Uj(*k - E j )^p 2m +2(X k ) 

d t X k = ±i =f 7T rr , (34) 

so if Afc = Ei at t = 0, it stays there forever, and it is consistent with the equations of motion to 
freeze some A& at the double roots of Q2n+2(A). This however cannot be done in an arbitrary way 
as we now explain. 

Let us assume that Q2n+2(A) has a real root at A = E. This means that ^4 2 (A) + -B(A)C(A) 
vanishes when X = E. But for real A, ^4(A) is real and one has C(A) = B(X) so that A(X), B(X), 
C(X) must all vanish at A = E. In particular, recalling eq, (|18p . this means that one of the separated 
variables, say Ai, is frozen at the value E, provided bb ^ 0. This implies also that X — E divides 
simultaneously A(X), B(X) and C(A), and therefore (A — E) 2 divides Q2n+2(A). A single real root 
of Q2n+2(A) is necessarily a double root. Writing: A(X) = (A - E)A(X), B(X) = (A - E)B(X), 
C(A) = (A - E)C(X), Q 2n +2(X) = (A - E) 2 Q 2n {X), we see that: 



Q 2 n(A) = H(X - e,) 2 U(X) 2 + B(X)C(X) 
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so we can apply the same discussion to the real roots of Q<2 n (ty- We deduce from this that the 
multiplicity p of the real root E of Q2n+2{X) is even and that (A — E) p ' 2 divides simultaneously 
A(X), B(X) and C(A). In this case, p/2 separated variables Ai,...,A p / 2 ar e frozen at the real value 
E. 

Nothing as simple holds in the case of a simple complex root. So, let us assume that there is 
a double complex root E. Of course the complex conjugate E is also a double root. What can be 
said in general is that the values of Ai,...,A n freeze by complex conjugated pairs. 

In fact we have A(A) = A 2 (X) + B(X)C(X) hence 

A'(A) = 2A(X)A'(X) + B'(X)C(X) + B(X)C(X) 

On a degenerate spectral curve, A(A) has at least one double zero E. Therefore 

= A 2 {E) + B{E)C{E) 

= 2A(E)A'(E) + B'(E)C(E) + B{E)C'(E) 

By contrast to the real case, we cannot infer from these equations that C(E) = 0. This is related 
to the fact that for critical spectral curves associated to unstable critical points, the corresponding 
level set of the conserved Hamiltonians has dimension equal to the number of unstable modes. 
So we expect that some separated variables A,- can remain unfrozen. But if one of them freezes 
at E, that is if C(E) = 0, the first equation implies A(E) = and the second equation implies 
B(E)C'(E) = 0. Therefore if the zero of C(X) at A = E is simple, then necessarily B(E) = 0. But 
B{E) is the complex conjugate of C(E) which must therefore vanish provided bb ^ 0. From this we 
conclude that another separated variable A& freezes at E. 

In contrast to the real case, freezing is not compulsory, but it is the possibility to freeze the Aj 
by complex conjugated pairs that leads to the description of the real slice and the stratification of 
the level set. 

3 The one-spin model. 

The model with a single spin coupled to the oscillator is interesting, because it illustrates most of 
the points discussed so far. In this case, there are only two conserved quantities, Hi and H2, which 
read: 

Hi = 2eiaf + bsf + 6a J" (35) 

H 2 = bb + sf (36) 

The rank of the momentum map can be either 0, 1, or 2. The later case corresponds to generic 
Arnold-Liouville tori of dimension 2. Let us now recall briefly the discussion of critical values of 
the moment map given in a previous work |llj . 
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3.1 Rank zero 

The rank of the moment map vanishes on the critical points given by b = b = 0, s 1 =0. Hence we 
have two points: 

s\ = eis, ei = ±1 

The corresponding values P = (H\,H2) are: 

P 1 (t) = [2e 1 a,a], P 2 (i) = [-2e lS ,-s] 

To determine the type of the singularities, we look at the classical Bethe equations (see section [9]) 

which read: 

2E + S€l = <£> 2E 2 - 2e ± E + se\ = 
E — ei 

The discriminant of this equation is e\ — 2se\. When the spin is down (point P2{\)), the discriminant 
is positive, the two classical Bethe roots are real and this is an elliptic singularity, in agreement 
with the general analysis of section [9l When the spin is up (point -Pi(t)) we have real roots when 
e\ > 2s (i.e. the singularity is elliptic in this case), and a pair of complex conjugate roots E, E 
when e\ < 2s (i.e. the singularity is focus- focus in that case). The image of the moment map is 
shown on Fig. [2] in the unstable case. We see that the stable critical point -Pi(t) is located at the 
tip of the image of the moment map, whereas the unstable one P2GI-) nes i n the interior of this 
domain. 

In the one spin case, eq. (fT3j) reads : 

Qa{X) -4A 2 + 4F 2 + 2-^ + -^ (37) 



(A- £l ) 2 A-ei (\-etf 

For the values Hi = 2e\e\s,Pi2 = e i s it is completely degenerate and the spectral polynomial has 
two double roots E\, E2 identical to the classical Bethe roots, and reads: 



Q 4 (A) = 4(A 2 - eiA + -e x ) 2 = 4(A - ^) 2 (A - E 2 f 



where 

se\ 
Ei + E2 = ei, E1E2 = —— 

Let us now describe in more detail the real slices corresponding to the two critical values of 
(Hi,H 2 ). Recall that 

af = (e 1 -Ai)(2Ai-Mi) (38) 

We express first the separated variable /ii as a function of Ai: 

(M -ei) 
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Figure 2: The image of the moment map in the case of one spin, (ei = —0.707), with one unstable 
critical point. The green point is the stable (elliptic) point P^i-l)- The red point is the unstable 
(focus- focus) point -Pi(t)- It is in the interior of the image of the moment map. The two boundaries 
of this domain correspond to spectral polynomials whith one double real zero, so that the moment 
map has rank 1 on the preimage of these curves. 



We have to distinguish the stable and unstable cases. In the stable case E\ and E 2 are real and 
we know that Ai has to be frozen at one of them, say E\. Then [i\ = 0, there is no problem of sign 
in eq. ([39|) and 

s\ = 2E 1 (e 1 - Ei) = 2EiE 2 = se 1 

as it should be. The level set in this case is composed of only one stratum consisting of the critical 
point itself where the rank of the moment map is zero. 

Let us now assume that we are in the unstable case, that is at the point i^Cf)) ei = +1 and 
2s — e\ > 0. In this case E\ = E, E 2 = E are complex and Ai remains unfrozen. Let us choose the 
+ sign in eq. (|39p . Inserting into eq. (|38p we obtain 



se\ 



so that, in the spin variables, we are precisely at the critical point. On this stratum of the level 
set the rank of the moment map is zero. It is very interesting to remark that the variable Ai 
completely disappears in this case (hence the rank zero). In terms of the separated variable, this 
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stratum of the level set seems therefore to consists of the whole Ai-plane. This is due to the fact 
that at the critical point bb = and the separated variables are not well defined and appear as 
singular coordinates. However, as we will show below, the Ai coordinate has the magic property of 
realizing a blowup of the singularity. 

Choosing now the — sign in eq. (|39j) yields: 



s{ = -sei -4Ai(Ai -ei) 

Writing that s\ is real, s\ = s\, gives 

Ai + Ai = ei (40) 

Setting Ai = x + iy, Ai = x — iy, we see that the real slice is given by the vertical line located at 
x = ei/2. It is not the whole line however. We find from bb = H2 — sf 

bb = 2s-e\- Ay 2 > 

The real slice of the reduced model corresponds in this case to the segment parallel to the imaginary 
axis delimited by the two double roots E and E of the spectral polynomial. This manifold has real 
dimension one. 

As we have discussed before |11| 118] . this corresponds to a real slice of the model reduced by 
the global U(l) action eq. (|2ip . From the viewpoint of the original model, we have to reintroduce 
this global angle and the real slice becomes a two dimensional torus pinched at the unstable critical 
point. 

Hence in this simple case we see that the level set of the unstable critical point is composed 
of two strata : the critical point itself where the rank of the moment map is zero and the pinched 
torus where the rank is two. 

The dynamics of the model on this torus is the composition of a large motion in which the 
oscillator amplitude bb goes to zero at time ±00 but reaches a finite maximum at a finite time, and 
a global U(l) rotation. Only the former movement is captured by the reduced model, and it maps 
into the finite segment of the Ai variable which we have just described. To see this, let us consider 
the equation of motion for the flow generated by H± : 



whose solution is: 



d tl \ 1 = 2i(\ 1 -E)(\ l -E) (41) 



Ai = ^^, X = X e^-^ (42) 

1 — A 



The reality condition eq. (|40p becomes: 

(E - E){X - X) = 

which imposes Xq = Xq so that Xq is real. Its absolute value can be absorbed in the origin of time 
t\. Only its sign matters. The constraint sf < s is equivalent to Xq < 0, and we recover the fact 
that Ai runs along the line interval joining E and E. 
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3.2 Rank one 

Physical configurations at which the rank of the moment map is equal to one correspond to a 
spectral polynomial with one double zero, which is necessarily real. We therefore write: 

Q 4 (A) = (A + ao) 2 (6 2 A 2 + 6iA + 6 ) (43) 

where a$ and bi are real. We denote: 

1 
a = - x - ei 

Taking into account the usual constraints on Q^X) saying that it depends only on two free pa- 
rameters H\ and H2 , see eq. ( [37]) , we can determine bo, b\, and 62 in terms of x by solving linear 
equations. We find: 



-4x, 


b = 


£l x 3 - 


efx 2 + 


s 2 




X 2 




i 
x — 


2eix 3 


-4s 2 








2x 






3x 4 - 


- 8eix 3 


+ 4ef x 2 


-As 2 





Hi = 

H2 = Ax 2 

These are parametric equations for a curve which coincides with the boundaries of the image of the 
moment map in the (-Hi, H2) plane, for which an illustration can be found on Fig. [2]. 

The determination of the real slice is straightforward in this case. Indeed, because we have a 
real double root, the separated variable is frozen at the value Ai = — ao, fJ-i = 0. This single point of 
the reduced model corresponds to the one-dimensional orbits under the global rotations generated 
by H2 characterized by a common phase on b and sf : 

s{ = --x(x-2ei) 



bb 



(2s - x 2 + 2eix)(2s + x 2 - 2e x x) 

Ax 2 
xb 



Hence, we have shown that the preimage of a point on the boundary consists of only one stratum 
which is a circle S 1 . Notice that when 

x 2 - 2eix + 2es = 0, e = ±1 (44) 

we have s\ = es and bb = and this corresponds to the points -Pi(|) and i^Ct)- 

In fact there is a simple relation between rank and rank 1. The spectral curves eq, (|43j) 
degenerate when the polynomial &2A 2 + b\\ + bo has a double real root. Its discriminant is: 

b\ - Ab b 2 = 16 ^ ~ 2 ^ X ~ 2s) f ~ 2 ^ X + 2S) = -64 bb < 

x z 
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It vanishes precisely when eqs. (|44p are satisfied. These equations are nothing but the classical Bethe 
equations eq. (|109p as can be seen by setting x = 2E. Hence in the stable case, the discriminant 
vanishes for real values of x, and the critical points are on the boundary lines of rank one. In 
the unstable case however the discriminant vanishes for complex values of x and the critical point 
cannot lie on the boundary. 

3.3 Reality conditions on a generic torus. 

In this simple case of a single spin it turns out that one can work out explicitly the exact real- 
ity conditions, for an arbitrary choice of the conserved quantities H\ and H 2 . For this, express 
everything in terms of Ai and \\. First, we have 

s+ = 26(ei - Ai), si = 26(ei - A x ), 

In the reduced model, H 2 is a constant and we have 

(H 2 -bb) 2 = (s\f = s 2 -stsl 

which yields 

{bbf + {AX - 2H 2 )bb + H 2 - s 2 = 

where X = (Ai — £i)(Ai — ei). Finally eqs. (j35|) and ([36]) give: 

H 1 = 2e 1 H 2 -2bb(\ 1 + \ 1 -e), (45) 

Eliminating bb between these two equations, we get : 

R : (F 1 -2e 1 i7 2 ) 2 + 4(^ 1 -2eiif 2 )(Ai + Ai-e 1 )(i? 2 -2AiAi + 2ei(A 1 + A 1 -ei)) 

+4(tf 2 2 -s 2 )(A 1 + A 1 -ei) 2 = (46) 

To obtain the equation of the real slice we ask that Ai and Ai are mutually complex conjugates. 
Setting Ai = x + iy, Ai = x — iy, the expression R is quadratic in y and cubic in x. It can be 
uniformized by Weierstrass functions. Setting: 

-6eiu + We 1 H 2 - 3#i - 2e\ -3v 



y 



4(-3n + 2H 2 - ef ) ' 4(-3« + 2H 2 - e 2 ) 

the reality condition for Ai takes the standard form for a planar cubic: 



with 



v 2 = 4u 3 - g 2 u - g 3 



g 2 = ^(H 2 + 2e 2 if 2 - 3eifTi + e\ + 3s 2 ) 
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g 3 = —(8H$ + 24eJ-H| - 36eiHiH 2 - 72s z H 2 + 24e?# 2 + 27 H( - 36e?#i + 36ef s z + 8e$) 
The range of physically acceptable values of x is determined by imposing: 

II ^i ~~ 2eii?2 . „ z rT 7, tt Hi — 2e\H2 , , 

tt = - 2(Al + Al - tl) ao,or Sl = g 2 -^ = g2 + 2(Ai + Ai _ ti) e{ - 8 . +s) 

where we have used eq. ([43]) . 

To help visualize these real slices in the Ai complex plane, we consider a small circle around 
the image of the unstable point in the (H\, H2) plane: 



Hi = 2se\ + r sin 9, H 2 = s + r cos ( 



(47) 



and we obtain a family of real slices shown in Fig. [3]. When Hi = 2eiH 2 , H 2 7^ ±s, the reality 
condition eq. (j4*6]) simplifies to 

Ai + Ai - ei = (48) 

In fact, on the line Hi = 2eiH 2 , the polynomial (54(A) depends on the variable A(A — ei) only which 
is invariant under the involution A — > ei — A. The slices thus become the vertical real line x = ei/2 
when 9 = 6 C where tan0 c = 2ei. 




Figure 3: The real slices when we run around the small circle around the singularity. €\ = —0.707, 
r = 0.3. 
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The two critical points lie on the line H\ = 2e\H2- When we are close to them, 

H\ = 2e\se\ + 5 Hi, H2 = se\ + a 5H\ 
Keeping the first order term in eq. (|46p we see that the real slice is the product of the straight line: 

Ai + Ai - ei = 
by the circle: 

(1 - 2eia)[2(Ai - ei)(Ai - e x ) - sex] - 2sei a (Ai + Ai - e) = (49) 

We can rewrite the above equation as: 

— Oi S61 

(Ai — A Q )(Ai — A a ) = [E\ — X a )(E2 — X a ), X a = ei + 



1 - 2ei« 
where Ei are the solutions of the classical Bethe equation: 

2E{E-ei) + se 1 = 

In the stable case e\ = —1, the Ei are real and when a varies, we have a pencil of circles with 
limit points located at X a = E\ and X a = E2, see Fig. [I]. 

In the unstable case obtained when e\ = +1 and e\ < 2s, the Ei are complex conjugate to each 
other and we have a pencil of circles with base points located at X a = E and A Q = E, see Fig. [5]. 

This expansion of eq. (06]) gives a very reasonable approximation of the real slice for (H\,H<2) 
close to the image of the unstable critical point -Pi(t) under the moment map. In Fig. [6], we see 
that for such values, the exact real slice is composed of two parts. One part is close to the vertical 
segment joining the two double roots E and E of the spectral polynomial associated to -Pi(t)- The 
second part is close to a circle belonging to the pencil just described. Here, we see a rather striking 
consequence of the fact that separated variables are singular in the vicinity of the critical points. 
Whereas the pinched torus containing the unstable critical point Pi (7) appears as a segment of the 
vertical line eq. (|48p in the Ai plane, any generic invariant torus arbitrary close to it is projected 
into a closed curve, whose shape depends crucially on the direction in the (H\,H2) plane along 
which we are approaching the critical value. The additional arc, well approximated by an element 
of the pencil of circles, corresponds to physical configurations which remain close to the unstable 
point -Pi(t)- To establish this, it is useful to use the quadratic normal form of the moment map in 
the vicinity of P\ (f) [H] , whose construction is recalled in section [9l 

3.4 Normal coordinates. 

To study the vicinity of a critical point, we can use the normal coordinates defined in eq. (1115p . We 
have 

«2 / V a l a 2 
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C(X) 



A-ei 



1.5 



0.5 



-0.5 



-1.5 




e-,/2 



Figure 4: The pencil of circles in the stable case e\ = —0.707 



The zero of C(A) is the separated variable Ai. We have 



A, 



E 2 o! 2 C x + E x a! x C 2 
a' 2 Ci + a[C 2 



(50) 



It is important to remark that since Ai depends only on the ratio C1/C2 a small neighborhood of 
the critical point is mapped to the whole Ai plane. 

We can parametrize the normal coordinates in terms of action-angle coordinates making explicit 
the reality conditions. 

In the stable case, we set: 

C j = A j e ie \ Bj = A je - ie J 

with real parameters Aj, 9j. The real slice in the Ai variables is therefore given by: 

E 2 a' 2 A 1 + E 1 a' 1 A 2 e iB 
Xl{d) - a' 2 A x + a[A 2 e^ ' e - e *-^ 

Taking the complex conjugate, we can eliminate between Ai and Ai. We find: 

(oi) 2 Al(Ai - £?!)(Ai - E{) - (4) 2 A?(Ai - E 2 )(\i - E 2 ) = 
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1.5 



0.5 



-0.5 



o c -t- jt/4 
6 C + 2ji/4 
6 C + 3ji/4 
6 C + 3.5jt/4 
6 C + 4.5ji/4 
6 C + 5jt/4 
6 C + 6jt/4 
6 C + 7ji/4 



-1.5 




Figure 5: The pencil of circles in the stable case (ei = —0.707. 



and we recover our circle pencil with limit points E\ and E 2 corresponding to A\ = and A 2 = 
respectively. 

In the unstable case, we set: 



B L = A lP e ie , C 1 = Aip~ 1 e- ie , 



Bo 



A 1P -V 9 , 



C 2 = A lP e 



-iO 



where p is a real positive number. Then: 



Mp) 



E 2 a' 2 A 1 + Exd^Axp 2 



a' 2 A 1 + a[Aip 2 
Taking the complex conjugate, remembering that E\ = E 2 , we can eliminate p 2 and get: 



I \2 All 



(ai)M^(Ai - £i)(Ai - Ek) - (a' 2 )M^(Ai - E 2 )(X l - E 2 



This is a circle pencil with base points E\ and E 2 . 
The conserved quantities are given by: 



H 2 



A 2 A 2 



,A\ 



-r + 



sA 2 



2(ei - £?i)oi 2(ei - £ 2 )a' 2 
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• / 

£1 

2 



Figure 6: The exact real slice (in blue) compared to the corresponding circle (in orange) of the 
pencil in the unstable case. The green dots are the branch points of the spectral curve. The 
vertical line is the line Ai + Ai = e±. The orange points are the classical Bethe roots, (ei = —0.707, 
#i = 2ei# 2 + 0.1, #2 = 1.) 



so that using the parametrization eq. (|47p . the equation of the pencil becomes: 

sin (9 [2(Ai -ei)(Ai-ei) -s]-2cos6> [2ei(Ai - ei)(Ai - ei) + s(Ai + Ai - 2e x )] =0 



It coincides with the pencil eq. (!49p as can be seen by setting a = cot 9. It is shown in Fig. [5]. 
Again, we stress that a specific circle of the pencil is specified by the direction in the (Hi, Hz) 
plane through which we reach the critical point. This is a typical blowup of a singularity. 

As we see on Fig. [6], if we consider a generic torus close to the pinched torus containing the 
unstable point Pi(t), the corresponding exact real slice in the Ai plane comprises a large circle and 
a segment very close to the line Ai + Ai — e\ = 0. This segment is absent from what we obtained 
from the quadratic normal form, which give us only the arc of circle belonging to the pencil. 

To understand this, we first emphasize that replacing the moment map by its quadratic normal 
form has a dramatic effect on the real slice corresponding to the pinched torus: instead of the 
segment joining E and E in the Ai plane, only the two end-points E and E survive after making 
this quadratic approximation. Indeed, in the immediate vicinity of the critical point, the pinched 
torus has the shape of two cones that meet precisely at the critical point. One of these cones is 
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associated to the unstable small perturbations, namely those which are exponentially amplified, 
whereas the other cone corresponds to perturbations which are exponentially attenuated. Taking 
into account the reality conditions, we can parametrize as follows: 

C 1 = A 1 , B 1 = A 2 , C 2 = A 2 , B 2 = A 1 

The action variables are: 

B X C X = AxA 2 , B 2 C 2 = A X A 2 

When we are sitting right on the pinched torus, these two action variables should vanish. We have 
two ways to do it corresponding to the stable and unstable cones: 

A x =0, A 2 + 0, or ii^O, ^2=0 

Inserting into eq. (|50p we see that they correspond to Ai = E\ and Ai = E 2 respectively. Hence the 
variable Ai remains trapped at these points in this approximation. 

Now, consider a solution of the equations of motion, and suppose we pick an initial condition 
close to the critical point Pi(t)j but not exactly on the pinched torus. Let us further assume 
that this initial condition is much closer to the stable cone than to the unstable one. The point 
representing the system in configuration space will first move towards Pi (j) along the stable cone, 
but after a finite time, while remaining close to Pi(t)> it will be ejected away from it along the 
unstable cone. This part of the movement takes place in the vicinity of -Pi(t) an d is well captured 
by the quadratic normal form. It corresponds to the arc of circle that connects the complex root 
E s associated to the stable cone to its complex conjugate E u associated to the unstable cone. This 
part of the trajectory has no counterpart when we are exactly on the pinched torus. 

After this, the point representing the system goes far away from Pi (f ) along the unstable cone, 
and the quadratic normal form is no longer valid. But of course, the two cones are connected and 
any initial condition located on the unstable cone gives rise to a trajectory which reaches eventually 
the stable cone in a finite time. This is manifested by the part of the trajectory in the Ai plane 
which remains close to the finite segment of the line eq. (!48p that connects E u to E s . This part is 
fully non-perturbative, in the sense that it cannot be captured by the quadratic normal form, but 
it is very well approximated by the large motion along the pinched torus. 

In a sense the separated variable Ai and the normal coordinate d play complementary role and 
are each well adapted to describe a different part of the trajectory. The need to glue these two 
qualitatively different parts on a trajectory close to a prinched torus is an important feature which 
played a crucial role in our semi-classical treatment of the model [18] . 



4 Critical torii associated to equilibrium points for general n. 

In this section, we study the level set of a critical point where the rank of the moment map is 
zero. In Fig. [I] these are the pre-images of the green (stable) and red (unstable) points. Hence, we 
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consider the critical points of the moment map (rank zero) . 

b = b = 0, Sj = sej, s • = 0, ej = ±1 
We know that they correspond to maximally degenerate spectral curves where: 

71 + 1 

Q 2 n+2(A)=4j](A-^) 2 (51) 

1=1 

From eq. (!33p we also know that the zeroes of Q2n+2{X) are the roots of the classical Bethe equation: 

n 
3=1 ° 

Notice that we have the identities: 



o Uk( e 3 ~ E k) _ ep v^ P _ , s 



The variables Hi are given by 



nj(Ai-£;) 

As we have seen in the one spin case, the choice of sign here plays a crucial role in the description 
of the various strata of the level set. 

We will study the strata of the level set by constructing the solutions of the equations of motion 
with generic initial conditions on the level set. The equations of motion eq. (120j) become 



We wish to solve these equations with the proper choice of signs and initial conditions so that the 
spins variables are real. 

We first reconstruct the spin variables. The s^ are given by: 

*;(*)= n 2 T rf ^'fl (53) 

iifc^( e j- e fc) 

where we have introduced the polynomial: 

V(X,t)=l[(X-X j (t)) 
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To reconstruct s|, we have to build the polynomial P n +i(A). From eqs. (l22|5ip . we have: 
P n+1 (\i) = ±VQ2n + 2(\i) = ±2 HiXi - Ei), P n+l (X) = 2(A n+1 - ai(e)A n + • • •) 

The choice of signs here is crucial. It is the same signs which appear in the equations of motion. 

If we take the + sign for all i, then obviously P n+ i(X) = 2]^(A — E{), and taking the residue 
at A = €j in P n+ i(X)/ n«(A ~~ e «)> we fi n d s j = se i- This is the static solution corresponding to the 
critical point. 

To go beyond this trivial solution, we divide the A, into three sets Xf G £ + , X® € £° and 
XJ G £~ depending on the sign in this formula. In particular £ is the set of Aj frozen at some E\. 
The equations of motion become: 

dtX r = - 2 M X * ~ El) , d t Xf = 2 l U ' l / (X i~ El) (54) 

where the prime means that elements in the set £° are excluded because they cancel between 
numerator and denominator. Therefore, for each Ei ^ £° we can write: 

y £ -y }* =-2,/ dz.u'^-E k ) = • a[{E) + El) 



XT A 



+ "i 



where C^ is a big circle at infinity surrounding all the Xj and we have defined T,^ = ^2jXj, 
°i(E)=Y!k E k- Hence: 

log ^|^y = -2iE t t -2if dtfo - a[(E)) 
where we have defined the polynomials: 

V ± {X)= n (A-AP), P (A)= n(A-A°) 

Remembering that 

b = -2ib(Y 1 ' l -a' l {E)-uj/2) 

we arrive at: 

V-(E l ) = b(t)X l V + (E l ), X l = X l (0)e- i ^+^ t , E t $ £° (55) 

Denoting h±,uq the number of elements in £^ \£° respectively, we get a set of n — no + 1 
linear equations for the n + + ra„ + 1 = ra — uq + 1 unknown coefficients of the polynomials V- (A) 
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and b(t)P+(X). As discussed in subsection 12.31 we expect that freezing uq separated variables Xj 
yields a stratum of dimension n — uq + 1 on the pre-image of the moment map corresponding to 
the degenerate spectral curve. 

Simple examination of the quadratic normal form of section [9] shows that the complex numbers 
2E[ + uj are the eigenfrequencies for the linearized equations of motion in the vicinity of the critical 
point. Eq. (|55p can therefore be interpreted as a non linear superposition of normal modes yielding 
a global motion along the critical torus. 

The next step is to build the polynomial P n+ i(A), using the fact that P n+ i(Aj) = ±2j^(Aj — E{): 

i x- Etfeo 11a^a-<A a i ) ' +i\ ) 

To derive this formula, we used the fact that the polynomials P n +i(A) and 2j^(A — E{) have the 
same terms of degrees n + 1 and n; the last statement comes from the relation £\. €j = ^ E\ which 
is a direct consequence of the classical Bethe equation. Similarly, we can also write: 

+4(A + S' x - a[(E))P + (X)P-(X)P (X) 

Note that the last term in the right hand side is necessary to adjust the terms of degree n + 1 and 
n in A between the two sides of the equation. Again, we used the relation Y2 j e j = Yli El- These 
two expressions for P n+ i(A) motivate the following definitions of polynomials 5±(A): 

s + w - £n<\--«fe^? 

so that we can write: 

P n+1 (X) = 2Y[(X-Ei)-4S + (X)P+(X)r Q (X) (56) 

i 

P n+l (X) = -2Y[(X-E l )+4S-(X)P-(\)Vo{X) (57) 

i 

Note that <S+ (A) has degree n_ — 1 and 5_(A) has degree n+ + 1. Now, we have: 

4 JJ(A - E t f - P n 2 +1 (A) = 16 S^(X)S + (X)P (X)P(X) = AbbP{X)P(X) 
i 
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where V{\) = Uj( x ~ x j) = P-WPoWP+W and V{\) = W^X - Xj) is the complex conjugate 
of V(X), that is P(\) = V(X). Therefore: 

bb P(X) = 4 5_(A)5+(A)P (A) (58) 

So the zeroes Aj of P(\) split into the zeroes A+, A~ and A° of 5+(A), 5_(A) and Vq(X) respectively. 
A direct consequence of these definitions and of eqs. (f56|) and (f57|) is that: 



P n+ i(Xf) = +2 J](A+-^) (59) 

i 

Pn+iCK) = -2]JCK-Ei) (60) 

I 

By the definition of the A^'s, we see that the set £° is self conjugate (see section [2T3|) and that 
"Po(A) = ^o(A)- The above definition of 5_(A) shows that the coefficient of its term of highest 
degree is equal to one (this is not the case for 5 + (A)). Because of this: 

S-(X)=P-(X) (61) 

Combining this with eq. (|58j) we get also: 

S + (\) = ±bbP + (\) (62) 

Comparing the terms of highest degrees in 5 + (A) and P+{X) gives: 

bb = AY U ff°[ K ~ E l ) (63) 

It is expressed only in terms of A~'s. So, if ra_ = 0, we recover the fact already mentioned that 
bb = and the system remains at the critical point. 

At this stage, we are ready to enforce the reality condition. As discussed in subsection I2.2( the 
real slice is obtained by imposing that the set {Aj} be the same as the set {A*}, where in the rest 
of this section we denote by z* the complex conjugate of z. Equivalently: 

p(\*) =V{X)* 

for any A. From the discussion in subsection 12.31 we know that the frozen variables A^ appear 
in complex conjugate pairs so that "Po(A*) = Vo(X)*. We also know that "Po(A) = Po(X) so that 
"Po(A*) = "Po(A)*. The above reality condition becomes then: 

P^{\*)P + {\*) = V-{\)*V + {\)* 
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It is clearly sufficient to impose simultaneously: 

p_(\*)=V_(X)\ P + (X*) = V+{\)* (64) 

But we now claim that this condition is also necessary. This comes from the fact already noted 
that the sign of P n+ i(X)/ \\i{X — E{) is positive for X = Xf 01 X = Xf and negative for A = A~ or 
A = \T . So the roots Xf have to be complex conjugates of Xf and likewise, the roots A~ have to 
be complex conjugates of A~. An interesting and useful consequence of this is that we must have 
deg V± = deg S± = deg P±, which requires n + = re_ — 1. Since n + + no + n_ = n we find: 

1, , 1, 

n+ = -{n-l-n ), n_ = -(n+l-n ) 

Because the coefficients of highest degrees of 'P-(A) and "P+(A) are set equal to one, the above 
constraints (|64p are equivalent to (assuming of course that P + (X) and V-(X) are mutually prime 
and similarly for V+(X) and V-{X)). 

P^(X*) __ fV-(X)\* 



V + (X*) \V + (X) 

We also note that we should add the constraint b = b* . This plus the fact that these two polynomials 
involve a total of n + +n_ unknown coefficients, shows that it is necessary and sufficient to enforce 
the following conditions for the n + 1 — uq = n + + n_ + 1 roots E\ which don't belong to £ : 



bV + (E;) \VP+(E, 

As we have seen, the general solution of the Hamiltonian evolution on the critical torus, eq. (|55p 
implies: 

To evaluate the left-hand side of the conditions (|65p . we set A = E\ in eqs. (|56| |57|) . which gives: 

P n+ i(£ z ) = -4S+(Ei)V+(Ei)V (Ei), P n+ i(Ei) = AS^E^V^E^V^Et), E x $ £° 

and therefore: 

S+iEtYP+iEi) = -S-{Ei)V-{E{) (67) 

Since 7 3 +(£'/) 7^ 0, and remembering eq. (|55p this implies: 

SJM _ 1 

S+(E t ) bX x 
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But using eqs. (fBTj) and (|62|) . we get: 

w> + (*;?) 4x f 

where we define the index I to be such that Ej = Ef . Given eqs. (|66|) and (|68|) . the reality conditions 
(|65p and the condition b = b* are satisfied if and only if: 

X?X T =~ (69) 

Let us check that the above conditions also imply the fact that s| is real for any j. For this, we 
have to satisfy the relation P n+1 (A*) = P n+ i(A)*. But, as we have just seen, conditions (|B"9"j) imply 
that A^ = (\\ )*. From this, combined with eqs. (|59p . (|60p and the fact that the polynomial 
]^(A — Ei) has real coefficients we see that: 



,{ ± )\*\ — p . cxW 1 



P n+1 ((AfO*)=^+i(A} 

But since: 

P n+1 ((A°)*) = P n+1 (A°)* = 

we get: 

Pi+i(A*) = P n+ \{\i)* 

for all the n separated variables. And because P n +\{\) = 2A n+1 — 2ai(e)A n + ..., these n relations 
are sufficient to impose the fact that P n +i has real coefficients. 

Together with the solution fj55[> of the Hamiltonian evolution, the explicit form eq. (|69|) of the 
reality conditions characterize completely a class of solitons, which have been first constructed by 
Yuzbashyan |12j . and called by him normal solitons. Applications will be given in sections [6] and [7] 
for systems with two and three spins. As expected, time disappears from these conditions so that 
they reduce to constraints on the integration constants: 

x,(0)* x T (o) = -i 

This characterizes a stratum of dimension re - no + 1 on the real slice. 

5 Case of n double zeroes. 

Let us consider now critical values of the moment map whose pre-image contains configurations 
where the moment map has rank one. These are the points on the green and red lines in Fig.pQ. 
In this case, the spectral polynomial Q2n+2(A) can be written as: 

n 

Q 2n+2 (A) = 4 p 2 (A) J](A - Ei) 2 , P2 (X) = A 2 + b x X + b 



4 = 1 
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We have seen that p 2 (A) is positive for real A. Hence 

A = b\ - 46 < 

If n r is the number of real zeroes, we must have n = n r + 2m because the complex Ei must come 
in complex conjugated pairs. As was explained in subsection 12. 3j, we must satisfy the relations: 

b l = 2{<y l {E)-a l (e)) (70) 

2ai(£) = 2<7i(e)-]T 



and 

This has the same form as eq. (|31|) and plays the role of a compatibility relation between the real 
parameters b\,bo. We get therefore a one parameter family of such spectral curves. As we have 
seen before, in the case of one spin, this family corresponds to the boundary of the image of the 
moment map in the (Hi,H2) plane. 
We set now: 

» = 4 + ^(A-A-') 

so that P2(A) becomes a perfect square: 

IA 



and: 



where: 



so that: 



p 2 (A) = ^(A + A^) 2 



f,AA\\ 2n+2 
Q 2 n +2 (A) = 4 Xl-i (A + A' 1 ) 2 J] (A - A' 1 - 2A,) 2 



1=1 



A, 



h + 2Ej 



V0^(X) = 2 I !LL± J (A + A- 1 )n(A-A~ 1 -2^)^Q(A) 

The sign of the square root is changed by the transformation A — > — A -1 . This transformation 
leaves A invariant and corresponds therefore to the hyperelliptic involution on the spectral curve, 
which has been uniformized in the complex plane A. 
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The points (X k , ±/i&) correspond to the points (A k , Q(A k )) and (— AT/ , Q(—A^ ) = — Q(A k )) in 
the complex A plane. Among these 2n points only n of them correspond to poles of the eigenvector. 
We denote them by A k and the other ones are their image by the hyperelliptic involution A^ = 
— A^T . As we know, when we have a double real zero Ei, one of the Aj's is frozen at Ei. Some other 
variables Aj's can be frozen by complex conjugate pairs on complex double roots. Let us denote by 
£° the set of uq > n r values A$ corresponding to the frozen roots. We have uq — n r even and as a 
result, n — no is also even. 

The equations of motion eq. (|20p become: 

A k __ yf\A\ ni(A*-Afc 1 -2A) 



A * 2 nUCAjk-Aj^-Ai + Af 1 ) 

where the ' means that the terms corresponding to frozen roots are excluded because they cancel 
between numerator and denominator. From this, we deduce: 

^Al-2A J A k -l 2 J * 

Introducing the roots of the denominators: 



B j = A j + ^A] + l, B] = -B^ = Aj -^Aj + l 
the solution of the equations of motion is: 

[[ k {A k - BA 



Note that the eigenfrequencies y|A]\ A 2 - + 1 = 2y/p^Ej) have a slightly more complicated relation 
to the double roots Ej than in the case of a spectral polynomial with n + 1 double roots studied 
before. 

We introduce now the polynomials: 

7MA) = n'(A-A fc ), n(A) = lT(A-A£) 

k k 

satisfying the useful formulae: 

P+(A") = (-A)- n+n »X\_A k Vl(A), Vl(A^) = (A)- n+ "»n / (A fc )- 1 V+{A) (72) 

k k 

With these notations, the solution of the equation of motion becomes: 

V + {B 3 ) = X, V + (B]) = X,BJ n+ ^ JjAfc PUJBj) (73) 
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These are n — no equations for the normalised degree n — uq polynomial 7-+ (A) which is therefore 
completely determined. Note that unlike the case of a spectral polynomial with n + 1 double roots, 
we have access only to the evolution of separated variables (the variables b, b do not appear in these 
formulae). In other words, this gives the dynamics of the reduced system, and the collective motion 
associated to the global U(l) action has to be determined separately. 

The next step is to find the polynomial P n +i(A), which we express in the variable A. To simplify 
the notation, P n+ i(A(A)) is a rational fraction in A which will be abusively denoted by P n +i(A). 
Note that with this notation A n+1 P n+ i(A) is a polynomial of degree 2(n + 1) in A. We have the 
constraints: 

P n+ i(A fe ) = P n+ i(A v k ) = Q(A fc ) = -Q(Al) 



By analogy to eqs. (1561) . (|57j) . we can write: 

A n+1 P n+1 (A) = A" +1 Q(A)-2S+(A)P + (A)P (A)P£(A) (74) 

A" +1 P n+1 (A) = -A" +1 Q(A) + 25„(A)^(A)P (A)^(A) (75) 

The polynomials Vq{A)Vq{A) are introduced because both P n+ i(A) and Q(A) vanish at those 
points. The coefficients of degree 2n + 2 of A n+1 P n+ i(A) and A n+1 Q(A) are identical, and the same 
is true for the coefficients of degree 2n + 1, because of the relation (J7DJ). This leads to: 

deg 5+ = n — hq, deg S- = n — hq + 2 

These two polynomials are determined by the n — uq conditions: 



and additional conditions at infinity. For example: 

n+l 

S^(A) = 2[^±L\ A«-«o+2 + ... 




From P n+1 (A) = P n+ i(A"), and Q(A) = -Q(A"), we get: 

5_(A) = JjA fc A n - no+2 5+(A r ') (76) 

k 

Note that the quantity \\ k A& is given by the constant term in the polynomial P+(A), which is 
reconstructed from the relations (|73p . Eq. (|76p implies that <S_(A) ~ A 2 when A — > 0. Furthermore, 
the known coefficient of highest degree in <S_ (A) determines the coefficient of lowest degree in 5+ (A) . 
The n — riQ values of 5+(A^) are then sufficient to reconstruct the polynomial 5+ (A) because its 
degree is n — uq. 
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Proceeding as before, we have: 

(r -\ 2n 
^P) bb A 2 V(A)P(A) 

where: 

V{A) = V + (A)Vl(A)Vo(A)V%(A), P(A) = P + (A)Pl(A)P (A)P^(A) 

Therefore: 

(/] — r\ 2n 
^ft) bb A 2 P(A) = S„(A)S + {A)V Q (A)V2(A) 

So we conclude that the frozen zeroes of "P(A) are the same as those of 'P(A), that is: 

Po(A)Pj(A) = ?o(A)?J(A) 

We have then: 

(/i — r\ 2n 
^Uj 66A 2 P + (A)^(A)=5_(A)5 + (A) (77) 

Hence, the non frozen zeroes of P{A) split into the zeroes A& of 5+ (A) and A^ of 5_(A) and these 
two sets of zeroes are related by the hyperelliptic involution due to eq.([76|). Notice that by eqs. (jT4"|) . 
(f75|) . we have: 

P n+1 (A fc ) = +Q(A k ), P n+ i(Al) = -Q(Al) 

which are indeed compatible. On the real slice, this will imply that A/% is the complex conjugate of 
A/;. Because 5_(A) is equal to A 2 times a polynomial of degree n — hq whose coermcient of highest 
degree is known, we get: 

(/TaT\ n+1 
^Uj A 2 ^(A) (78) 

Substituting this equality in eq. ([77|) . we find: 

/ /i — r\ n ~ * 
S+{A) = l\ 4 ~ bbP +( A )> ( 79 ) 



Inserting into eq. (|76p we find the compatibility relation: 



4 niA t n t A* 

where have used the fact that n — no is even so that (— l) n_n ° = 1. 
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Let us now turn to the real slice. Reasoning in the same way as in section 0J it is necessary and 
sufficient to impose the n — no conditions: 

P+(Bf) (P+(Bi) 



V + (B r >*) \P+(B?) 
which are equivalent to (recall that ^(O) = n Afc) : 

P+(Bf) _ ( V + {B l 



(80) 



From eq. (fT3"j) . we have: 

x _ m -n+n (gl) 



V+{Bt 



V + (0)VZ(Bf) 
To express the right-hand side in eq. (]80j) . we write eqs. (|74h . (175p for A = B[, where I £ £°: 

^" +1 P n+1 (^) = -2S + {B l )V+(B l )V (B l )V^(B l ) 
5f +1 P n+1 (B z ) = 2S-{B l )Vl{B l )Vo{B l )V%{B l ) 

and therefore: 

S + (B,)P+(B,) = -S_(B,)P!(B,), /^° (82) 

Remembering eq. (|73p this implies: 



5_(B,) = (.^n-no+l^-n+no^ JJ Afc S+ ^ 

k 

Using eqs.(|75|). ([79|) . we arrive at: 



(83) 



P+(0)P+(BO X, 

Taking into account eqs. (f8Tj) and (f83j) . the reality conditions ([80]) become: 

XfXi = -(Bf) 2n - 2no+2 (84) 

Time disappears from these conditions which reduce to constraints on the integration constants: 

Xi(0)* X0) = -[B*) 2n - 2no+2 , I $ £° 



Together with the solution (|73|) of the Hamiltonian flow, these equations characterize completely 
any stratum of dimension n — tlq + 1 on the real slice where uq separated variables are frozen 
on the double roots of the spectral polynomial. They correspond to the anomalous solitons in 
Yuzbashyan's terminology |12j . 
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6 The two-spins model. 

We now give some details on the two spins model. Let us first write explicitly the Hamiltonians: 

H x = 2eisf + bsf + bs^ + ' 



H 2 = 2e 2 s 2 + bs 2 + bs 2 



ei - £2 

Sl ■ S~2 



ei -e 2 

H 3 = bb + sl + s z 2 

where s{ ■ s*2 = s\ s 2 + 2 {si s 2~ + s i~ s 2~)- The singular points are given by b = 6 = 0, s 1 = s 2 = 
so that we have four of them: sf = ±s, s| = ^s- The corresponding values P = {H\,H 2 ,H^) 
are: 

PlitA) = [^ + 26 lS , 2e 2 s-^, 2a] 

P 2 (t,4) = [2€ia- i ^, ^-26 2S , 0] 

P 3 a,t) = [-2e lS -^, -^- + 2e 2S , 0] 

A(l,4) = [^-26 lS , -2e 2S -^, -2 8 ] 
In order to determine the type of the singularities we write the classical Bethe equations: 

a(E) = 2E + -^ + -^ = (85) 

ej — t\ tii — e 2 

These are polynomial equations of degree 3. The type of the singulariy is determined by the number 
of real roots of these equations. If we have three real roots we have a stable elliptic singularity and 
if we have one real root we have an unstable focus- focus singularity. We easily see that P4Q,, 4) is 
always stable. See [11] for a detailed discussion of the parameter space (ei, e 2 ). 

6.1 Image of the moment map. 

The spectral curve reads in this case: 



Q*W -41*4.477 4- 2Hl 4- 2H2 ■ s 1 s (m 

(A-ei) 2 (A -e 2 ) z A - ei A - e 2 (A-ei)^ (A - e 2 )^ 

The image of the moment map can be drawn from the degenerations of this curve. An example 
is shown in Fig. pp. In this example the parameters (ei,e 2 ) have been chosen such that the points 
P 2 (t, 4), Pi(4, 4) are stable while the points Pi(|,t) an( A P3(i,T) are unstable. 

The spin variables are reconstructed from the separated variables (Aj, /Xj), i = 1,2. In particular, 
we have: 

ofa- Ai)( Cl - A a ) (Ai+Aa _ ea) A2- ei Ai- ei (87) 



ei-e 2 (Ai - A 2 )(ei - e 2 ) (A 2 - Xi)(e x - e 2 ) 

Ai)(e 2 — A 2 ), , A 2 — e 2 Ai — e 2 

(Ai + A 2 - ei) - 1/1— — - - v 2 — — 

e 2 -ei (Ai - A 2 )(e 2 -ei) (A 2 - Ai)(e 2 - ei) 
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where Vi = (A; - ei)(Aj - e 2 )^. 

6.2 Rank zero. 

The critical points are the points where the rank of the moment map is zero. They correspond to 
the complete degeneracy of the spectral curve: 

f = a\X) 

where a(A) is given eq. (|109p . The pre-images of these points are reduced to a point, in the case of 
a stable critical point, or to a pinched torus in the case of an unstable critical point. In that case 
non trivial solutions of the equations of motion exist which we now study. 
So, let 

Q 6 (X)=4(X-E ) 2 (X-E) 2 (X-E) 2 



where Eq,E,E are the three roots of eq. (!85p . We assume that Eq is real E is complex and E is the 
complex conjugate of E. We have the following relations coming from eq. (|85p : 

Eq + E + E = ei + e 2 
EE + E E + E E = eie 2 + |(ei+e 2 ) 

s 
E EE = -(eie 2 + e 2 ei) 

Since there is a double real zero, we must freeze one point of the divisor : 

A 2 = E , fi 2 = 

then there remains only (Xi,fi\) with 

^ ^ - Ep)^ - E)(X 1 - E) 
(Ai -ei)(Ai -e 2 ) 

then from eqs. (|87l) . (J88j) we get immediately: 

si = -2 e -^^((X 1 + E -e 2 )(X 1 -e 1 )T(Xi-E)(X 1 -E)) 



ei - e 2 v ' / 

si = _ 2 e -l—^-((X l + E -e 1 )(X l -e 2 )T(X 1 -E)(X l -E)) 
e 2 - ei V J 

The upper sign leads to the constant solution s\ = sei as we already know. So we choose the lower 
sign. Reality of s\ requires: 

A1+A1 = E + E 
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and we can write: 



-2 £l E ° ( - 2AiAi + ei(e 2 - E ) + EE) 
ei - e 2 \ 

-2 e2 ~^ f - 2AiAi + e 2 (ei - £ ) + ££ 
e 2 — ei V 



The equation of motion becomes: 

dAi 



2i{X 1 -E){X 1 -E) 

at 



This is exactly the same equation as eq.(|4Tj). whose solution is given by eq.(|42j). It reads: 

Xl = E ~ E ^ , X = X(0)e 2 ^ E - E ^, X(0) reaKO. 
I — X 

These solutions exist on the pre-images of the red and blue dots in Fig. [1] . 

6.3 Rank one. 

The lines of rank one correspond to the following degeneracy of the spectral curve: 

Q 6 (A) = 4(A 2 + aiA + a ) 2 (A 2 + &iA + b ) (90) 

We have four coefficients and three conditions on Qq(X). Hence we have a dimension one manifold 
of solutions. The coefficients bj are completely determined and there is one constraint between 
(ao,oi). Note that the case of rank is obtained as a special case of rank 1, when the polynomial 
A 2 + &iA + bo has a doubly degenerate root, which is necessarily real. Let us parametrize (ao,oi) 
in terms of (x, y) as follows: 

ay + e 2 x - 2eie 2 y + x-2e 2 -2ei 
ao = 2 ' ° 1 = 2 

Imposing the vanishing of the coefficient of A in eq. (l86|) . we find: 

h = - (y + x) 
Imposing that the coefficient of the double pole at A = e 2 is s 2 we find: 

, (e 2 2/ 3 + e 2 xy 2 -e^y 2 + s 2 ) 
»0 = 5 

yZ 

Imposing that the coefficient of the double pole at A = ei is s 2 we find that the parameters x and 
y are tied together by the relation: 

Si : s 2 {y 2 - x 2 ) + (ei - e 2 )(x + y - e x - e 2 )x 2 y 2 = (91) 
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This is nothing but eq. (|7ip when we set: 



e\s e 2 s 

y 



Once the spectral curve is known, we can read the values of the Hamiltonians Hi by writing it 
in the form eg. (186 j) . They read: 



1 2s 2 

Hi = --x 3 - -yx 2 + e ± x 2 + - -(y + 2e 1 -2e 2 ) (92) 

2 2 (ex - e 2 )x 

1 1 2s 2 

H 2 = -^y 3 -^xy 2 + e 2 y 2 -- -(x-2 ei + 2e 2 ) (93) 

2 2 (ex- e 2 )y 

3 1 s 2 s 2 

Hi = --(x + y-ex - e 2 ) 2 - e 2 x - exy + eie 2 + -(ex + e 2 f + ^ + —^ ( 94 ) 



Together with eq. (|91|) these are the parametric equations of the lines of rank one of the moment 
map. 

As we have seen, the discriminant b\ — 46o&2 < 0. It is zero when we are at a critical point, in 
which case A 2 + &iA + bo = (A — Eq) 2 with Eq real. But as soon as we leave the critical point it 
becomes strictly negative so that: 

X 2 + b 1 X + b = (X-E )(X-E ) 

with Eq and Eq being complex conjugate. The difference between the green lines and the red lines 
in Fig.[T] comes from the sign of the other discriminant a 2 — 4ciq. This leads to two very different 
situations which we describe below. 

First when a 2 — 4ao is positive the polynomial A 2 + a±X + ao = (A — Ex)(\ — E 2 ) has two real 
roots, so that Qe(X) has two double real roots. According to the general discussion we must freeze 
two separate variables on these roots, and hence Ai and A2 are completely frozen. This corresponds 
to a static solution of the reduced model and therefore to a Liouville torus reduced to a circle S* 1 
in the original model. This circle can be easily described. We can compute the components s 2 - of 
the spins by setting (Ai = Ex,/J,x = 0), (A 2 = £2,^2 = 0) in eqs.flU]), flEEl- We find: 

sf = 2 (ei ~ .^ l)(ei - E2) (Ex +E 2 - e 2 ) = \x(2ex -x-y) (95) 

(ex — e 2 ) l 

4 = 2 {e2 -^ ){e2 - E2 \ Ex + E 2 -ex) = \y(2e 2 -x-y) (96) 

[e 2 — ex) l 

The spins simply have a uniform rotation around the z-axis. The quantity bb is constant and can 
be calculated from bb = H 2 — s\ — s\. This is the situation that prevails when we move along a 
green line in Fig. pp. 
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The other situation occurs when a\ — 4ao is negative. The polynomial A 2 + a\\ + ao = (A — 
E){\ — E) has now two complex conjugated roots, so that Q%{X) has pair of complex conjugated 
double roots. This is the situation when we move on a red line in Fig.pQ until the discriminant 
eventually vanishes. At this point the two roots (E,E) become a double real root (so that Qq(X) 
has a quadruple real root there). Beyond that point we have two real roots and the color of the 
line turns green. 

The pre-image of a point on a rank one line can be conveniently understood in the vicinity of a 
critical point by using normal forms. Around a stable critical point, we have three elliptic normal 
modes. At the critical point the three action variables corresponding to these modes are equal to 
zero. When we leave the critical point on a rank one line, two action variables are kept equal to 
zero, and a third one becomes non zero. We get the three lines leaving the critical point from the 
three possible choices of the non zero action variable. The pre-image of a point on these lines is a 
circle. 

Around an unstable point, we have one focus-focus mode and one elliptic mode. The neighbor- 
hood of the critical point on the two dimensional face on which it lies is obtained by the activation 
of the action variables of the focus-focus mode, keeping the action variables of the elliptic mode 
equal to zero. On the other hand, the rank one line leaving the critical point corresponds to the 
action variable associated to the elliptic mode becoming non zero while the two action variables 
associated to the focus-focus mode remain zero. In phase space, the pre-image associated to the 
elliptic mode is a circle, while the pre-image associated to the focus-focus modes is a pinched torus. 
Hence the pre-image of a point on the rank one line is the product of a circle and a pinched torus. 

As long as a 2 — 4ao is negative, there exist non trivial solutions of the equations of motion. So, 
let us concentrate on that case: 

Q 6 (A) = 4(A - Ef(\ - Ef{\ - E )(\ - E ) 

One simple solution in this case is to freeze the separated variables on the pair of complex 
conjugated double roots (Ai = E,/j,i = 0), (A2 = E,fj,2 = 0). This static solution corresponds 
precisely to a circle in the original model, and eqs. (|95l96p are still valid in this case. This is the 
stratum n = 2, hq = 2 in the general theory. 

But we know that there exists another stratum, n = 2,hq = 0, of dimension 3 in the original 
model and 2 in the reduced model, corresponding to big motions on the pinched torus mentioned 
above. To uniformize the spectral curve, we set as before: 

A = ^^ + ^^(A - A"*) = -| + >/E(A - A-) (97) 

so that: 

(\-E)(\-E) = ^(A*-AA-l)(A*-AA-l) 



where we have set: 

A 



16 A 2 
61 + 2E - &! + IE 
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We introduce the coordinates (Ai,A2) associated to (Ai,A2) and the polynomial P+(A) = (A 
Ai)(A — A2). The solution of the equations of motion is given by: 



(Ai-fli)(A 2 -Bi) = Xi(*)(Ai-J3?)(A 2 -J3?) 
(A 1 -i*i)(A 2 -.B I ) = X [ (i)(A 1 -JB?)(A a -.B?) 



(98) 
(99) 



where: 



B x = A + y/A 2 + 1, B\ = A - \A4 2 + 1, P T = A + \/A 2 + 1, B\ = A- VA 2 + l 



and: 



Xl (t) = A 1 (0)e^^+ Tt , Xi(t) = X-^Oy^l"^ 11 
The reality condition reads: 



X 1 (0)X I (p) = -(B 1 ) 



(100) 



We can express the quantity bb : 



bb 



IAI 



1 



4 AiA 2 AiA 2 



The spins can be reconstructed through the formula sf = "+u e jJ 

F & 3 Uk^i( e J- e k) 



polynomial P n+ i(A). Adding eqs. ()74|75p and using eq, ([76|) and eg. ([72]) . we get: 



For this we need the 



A 3 P 3 (A) 



IAI 



A 3 



AiA 2 AiA^ 



[A 3 P + (A")P + (A") + (A") 3 P+(A)7> + (A)] 



hence: 



^s(A) 



1 
~PP 



1\2 



PP(A - A~y + (PS + PS) (A - A~ x ) 



where: 



+(SS + 3PP + P + P - 1)(A - A" 1 ) + 2(PS + PS + S + S) 



S = A! + A 2 , P = A!A 2 , S = A!+A 2 , P = A 1 A 2 



From eq. (|97p . we see that this is polynomial of degree 3 in A, as it should be. 

One can obtain the reality condition expressed directly on (Ai,A2). Eqs. (|98|99p are a linear 
system for the symmetric functions S, P. One can eliminate X and X between these equations and 
their complex conjugate, taking into account eq. (|100p . The result can be written in the form: 



s s nf . j, r n/ 1 \ s s A . j 

p + T = -2(A + A), l + _ l + ? + _ ?= 4AA 



101) 
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An example of the motion of the variables Ai, A2 is shown in Fig. |7|8| . and the associated motions 
of sf , s| and bb are shown in Figs. (j9ll0p . We recover the soliton-like behavior of the large motion 
already discussed for the rank zero case. One difference is that now, the spins are no longer along 
the z axis as times goes to ±00. But we also see that the pulse in the oscillator energy bb(t) develops 
a non-trivial time dependence, with oscillations superimposed to a typical solitonic enveloppe. This 
reflects an important qualitative difference between the rank zero and the rank one case. In the 
former case, the pre-image is just a two-dimensional pinched torus, whose small cycle corresponds 
to the global U(l) action generated by H3. This cycle is not seen in the reduced model, for which 
the large motion connecting the unstable manifold to the stable one is one-dimensional. But for a 
point on a line of rank one, the pre-image is the product of a circle by a two-dimensional pinched 
torus. As a result, the winding motion around the small cycle of this torus becomes observable for 
the reduced model, as manifested by Figs. (|9llOI) , 
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Figure 7: Motion of the point say Ai of the divisor. The second point A2 is hidden into the black 
dot. (ei = -1.2,e 2 = -1.735). 
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Figure 8: A zoom of the motion of the divisor Ai, A2. (ei = —1.2, €2 = —1.735). 




Figure 9: The components s\ and s\ of the spins, (ei = —1.2, €2 = —1.735). 
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Figure 10: The quantity bb. (a = -1.2, e 2 = -1.735). 
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7 The three spin model. 

7.1 Critical stratum of dimension 4 

Let us examine in some detail the example of three spins, for which we focus on the case of fully 
degenerate spectral curves associated to the critical points: 

b = b = 0, s ■ = 0, Sj = ejS, ej = ±1 

Therefore, the spectral polynomial can be written as: 

4 

Qs(X) = 4l[(X-E l ) 2 
i=i 

According to the general analysis presented in section [U the different strata in the associated 
degnerate torus are obtained when uq = and no = 2. Consider first the case no = 0, (or in the 
general case no = n — 3). We have: 

deg V + = deg S + = 1 
deg V- = deg 5_ = 2 



Wc set: 



Equations (l55j) become: 



^(A) = A 2 -(TiA + a 2 , P+(A) = A-A 



-t — Ei ^=- + E l j + Xi \q — EiXi, 
b b b 



or, in matrix form: 



/! 



£1 



£? 



\1 



-©2 ^2 

E4 E 4 



X 1 \ 

x 2 

X3 

X A J 



We introduce the matrices: 



Mn 



/I ^1 S? 

1 E% E 2 

1 ^3 Si 

\l E A El 



Mo 



(E 1 E\ X 1 E 1 X 1 \ 

-E/2 E 2 X2 E2X2 

E3 E 3 X3 E3X3 

\E4 E 4 X4 E4X4J 



X 1 \ 

x 2 
x 3 

X 4 J 



( °2/b \ 
■<ri/b 
1/6 
Ao 



V 



M, 



Mi 



1 
1 

Vi 



1 = 1---4. 



-©2^2 
£z3AL3 



^1 

E?2 

E3 

E 4 



El 
El 
El 
El 



(I E\ X 1 ExX^ 

1 E 2 X 2 E2X2 

1 E 3 X3 E3X3 

\ 1 E A X4 E4X4 J 



EiX ± \ 
E2X2 
E3X3 
E4X4J 



M A 



/l E x X x E ± XA 

1 E% X2 E2X2 

1 E3 X3 E3X3 

\ 1 E4 X4 E4X4 J 
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Denoting Di = det Mi we have: 



x Dl D * D * T, D ° nn>>\ 

Ao = ^' CT2 = A' ai = W b = ~Dl (1 ° 2) 



Moreover, eq. (|63j) becomes: 



- bb = 4 / Q(X-)-Q(\-) ^ + \~Q(\~) - X^Q(X 2 ) 



where we have introduced the polynomial Q(\) (not to be confused with Qs(A)) defined by: 

4 

g(A) = IJ(A--Bj) 

The right hand side in the expression of 66 is a function of a\ and 02 since AJ" and A^ are the roots 
ofP_(A). 

In Fig.|ll| we show an example of parameters e-,-, ej such that the classical Bethe equation has 
four complex solutions. Note that if we choose for our critical point the ground-state of the diagonal 
part of the Hamitonian V ■ 2ejSj, we have ej€j negative for all j, and a simple graphical construction 
shows that the Bethe equations have at most a single pair of complex conjugate roots. In this case, 
according to our general discussion, uq = n — 1 and the preimage is composed of a two dimensional 
pinched torus, independently of the value of n. In the example of Fig.jllj. (ei, e2, 63) = (1, —1, 1) 
and (£1,62,63) = (—3,-2.7,0.5) so that only e\e\ is negative, and the critical point is a doubly 
excited state of ^ ■ 2tjSj. 

From the general form of the time evolution (I55D . we see that: 

and the reality conditions eq. (|69|) imply then that: 

Xi(t) = -4x,(o) 2 *K-t) 

As a result, if Xf(0) = 1/4 then we have bb(—t) = bb(t) so that bb has an extremum at t = 0. 
Only the relative signs of AQ(0) matters. In Fig. [12] we draw the function bb{i) when -Xi(O) = 
X 2 (0) = 1/2. The maximum is at 66(0) = \E X + E 2 - E\ - E 2 \ 2 . The time evolution of bb{t) 
shows a rather rich internal structure on top of an overal solitonic shape. Note that this trajectory 
lies on a stratum of dimension 4, which can be described, at least in the vicinity of the unstable 
critical point (ei, 62,63) = (1, — 1, 1) where the system is equivalent to its quadratic normal form, 
as a product of two pinched torii of dimension 2. The complexity of the motion is illustrated in 
Fig.pS] which shows the trajectory Ao(i) in the case -X"i(0) = ^(0) = 1/2. We see that it runs 
from one of the roots Ei to its complex conjugate. Because the dimension of this stratum is larger 
than two, it is quite easy to change the shape of these solitonic pulses. An other example is shown 
in Fig. Q3] when Xi(0) = -X 2 (0) = 1/2. The extremum at 66(0) = \E 1 - E 2 - E x + E 2 \ 2 is now a 
local minimum. 
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7.2 Critical stratum of dimension TV 

These observations can be extended for a system with an arbitrary number of spins. Let us consider 
a stratum of dimension N containing a critical point, where N is the number of roots which don't 
correspond to any frozen separated variable Aj. We have N = n-\-l — uq = n_ + n + + 1 and N is 
even because n_ — n+ = 1. In general, b can be written as the ratio of two determinants of size N, 
that is b = (-l) n +D /D n+1 where: 

D = I f\ E A E 2 ■ ■ ■ f\ E n - A X A XE ■ ■ ■ A XE^ n +^ l) (103) 

D n+l = \AEAE 2 ---A E {n — l) AX AXE--- A XE n + (104) 

Notice that the degree in energy of these determinants are: 

n_(?i_ + 1) n + (n + -l) 
d ° = 2 + 2 

n^{n^ — 1) n + {n + + 1) 
d n +i = n 1 7, 



so that 

do - d n+ i = n_ - n + = 1 

At the symmetric point Xi(t = 0) = -Xj(O) = 1/2, we have: 

Dn+i=Y[(E i -E j )(E i -E j ) 

i<j 

Dq = D n+1 Y,{E i -E l ) = 2iD n+1 J2 ImEi 

i i 

All ImEi are positive, so assuming that they are of order one: 

66(0) > iV 2 |min ImEi\ 2 
and we recover the super radiance phenomenon of Dicke. 

7.3 Critical stratum of minimal dimension 

We consider now the case with no = 2 frozen roots. This gives a smaller stratum of dimension two. 
In the general case we would take no = n — 1. We have: 

deg V+ = deg S + = 
deg V- = deg <S_ = 1 
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Figure 11: The 3 spins model. Determination of the parameters ej such that the classical Bethe 
equation has 4 complex roots. The blue curve is y = se±/(x — e±) + se2/(x — €2) + se^/{x — £3), the 
red curve is the straight line y = —2x. (e\, e2, 63) = (1, — 1, 1), (ei, £2, £3) = (—3, —2.7,0.5). 



We choose to freeze (A2,As) at (E2, E2). So Ai is the only separated variable which is not frozen. 
We can write: 

P_(A) = (A-Ai), V+(X) = 1, V (X) = (X - E 2 )(X - E 2 ) 

The conditions eq. (l55l) give 

V-(E 1 ) = bX 1 , V-(E l ) = bX l 

which are solved by: 

E\X\ — E\X^ — E\ — E\ 



A, 



X l -X l 



Xi -X T 



where: 



Ai(0)li(0) 



1 



The last equality yields the familiar compatibility relation: 

Ai + Ai = E\ + E\ 



(105) 
(106) 



(107) 



Since the unfrozen pair can be (Ei,Ei), i = 1 • • • n, we have n such solutions. We show an example 
of the small strata in Fig.|15j. We recover a simple solitonic pulse, in agreement with the general 
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Figure 12: The 3 spins model. bb(t) as function of time. X\(0) = ^(0) = 1/2 



idea, discussed in subsection 12.31 that in this case, the system dynamics can be mapped into the 
one of an effective model with a single spin. 

7.4 Relation with Normal coordinates. 



A natural question is to compare the exact solution Eq. ()102j) to the normal modes expansion 
Eq. (1116p . Normal modes are defined near a critical point which is precisely the pinch of the torus. 
So, we have to examine the exact solution near the pinch, that is to say when time t — > ±oo. 

As we describe in section[9l normal coordinates are of the form Bi = B(Ei), C\ = C(Ei), where 
E{ are the n + 1 double roots of the spectral polynomial. There, we show that: 

d(t) = Ci{0)e i(u+2E ^ 

This implies that Ci(t) remains proportional to 1/Aj(t), and we can choose the initial condition so 
that 

Ci(t) 



1 



hence, using Eq. (|116|) : 



b(t) 



V 



Xi(t) 
1 1 



*?<*&) 
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Figure 13: The 3 spins model. The trajectory of Aq. 



When t —> oo, Xi(t) — > or oo. Since bb remains bounded, terms with Xi(t) — > cannot appear 
when we take the limit of an exact solution. So, when t — > — oo only terms with Im Ei < 
contribute, and when t — > +oo only terms with Im Ei > contribute. By definition, we call £"j 
those classical Bethe roots with Im Ei > 0. Then Ej = Ei has Im E^ < 0. 
The simplest case is given by the stratum of minimal dimension Eq. f|106|) : 



'mm 



(0 



Ei - E, 



Xi(t)-Xj(t) 



Then 



■Tnin 



(t)\i 



Xi(t) 



y mm 



(*)l 



t—>+oo 



E i - E I 
Xi(t) 



So, we find the remarkable result that ^ m j n (£) is an exact solution which interpolates between the 
single normal modes Cj and Cj at t = ±oo. In terms of separated variables, this corresponds to 
n— 1 variables Aj frozen on n — 1 of the Ei, and the last Aj moving between Ei and Ej. As t — > ±oo, 
Aj — )• Ei, Ej and the normal modes correspond to all Aj frozen. 

More generally, when we consider the limit at t — > — oo in Eqs. (|103ll04p . we should keep the 
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Figure 14: The 3 spins model. bb{t) as function of time. Xi(0) = 1/2,X 2 (0) = -1/2 
dominant terms, which amounts to replacing (assuming for simplicity that all roots are complex) 



X 







\*5+r/ 



Then: 



b(t)\ 



t—t—oo 



E 



X\^{Ej-E 3 )X-{ty 



b(t)\ 



t— >+oo 



E 



U-jJEj-Ej) 1 



Hence, the solution appears as a superposition of normal modes at t — > — oo and the above formulae 
determine the asymptotic behavior at t = +oo. Of course, in between, the solution is a non linear 
superposition of normal modes which is completely beyond the scope of the quadratic normal form, 
because it describes a motion that starts from the unstable manifold on the critical torus and 
eventually returns to the stable manifold, after a big motion on the torus. The quadratic normal 
form is on the other hand very useful to describe trajectories in the vicinity of the critical point 
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Figure 15: The 3 spins model. The curve bb(t) on the small stratum. ^(0) = 0.5. 



but on a different pre-image of the moment map, which begin close to the stable manifold, and get 
kicked towards the unstable one. As for the single spin case, a generic trajectory on a torus close 
to the critical one will consist in two parts, one well described by the quadratic normal form, and 
the other well approximated by the generalized solitonic pulses constructed in this section. 

The various strata of the real slice are characterized by the number of modes appearing in the 
expansion of the solution at say t = — oo. In principle this could be controlled by sending some 
initial conditions Xj(0) to oo. But because of the reality conditions, this is a rather singular limit, 
and that is why an independent description of each stratum is necessary. 

8 Conclusion 

In this work, we have analyzed in detail singular torii in the classical integrable Dicke-Jaynes- 
Cummings-Gaudin model which describes a system of n inequivalent spins coupled to a single 
harmonic oscillator. These singular torii correspond to critical values of the moment map con- 
structed from the conserved quantities under the Hamiltonian dynamics. The level sets associated 
to such critical values appear to have a natural stratification, the dimension of each stratum be- 
ing equal to the rank of the differential of the moment map. Although the complexified model is 
the natural stage for the algebro-geometric method used to solve Hamilton's equations of motion, 
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from a physical perspective, it is crucial to look at these level sets on the real slice, whose precise 
identification has been one the main goals of this work. 

The dimension of a stratum on a singular level set is given by n + 1 — hq where no is the 
number of separated variables [15] which are frozen on some double roots of the associated spectral 
polynomial Q2n+2(A). The values of no vary among the integers ranging between the number of real 
double roots of Q2n+2(A) and the total number of double roots of this polynomial, given the fact 
that frozen roots occur in complex conjugate pairs. The fact that the possible dimensions of strata 
belonging to a given level set vary by multiples of 2 has a simple interpretation. We know that 
singularities in the Dicke-Jaynes-Cummings-Gaudin model have a quadratic normal form which can 
be either of elliptic or focus-focus type |llj . The quadratic normal form gives a picture of a stratum 
on a level set that is, in the vicinity of a critical point, the product of an m-dimensional torus and p 
two-dimensional pinched torii. This gives a stratum of dimension m + 2p where p ranges from zero 
to the number of distinct focus-focus singularities appearing in the quadratic normal form, which 
is equal to the number of complex conjugate pairs of double roots of the spectral polynomial. The 
total number of double roots of this polynomial is then equal to n + 1 — m. 

We gave a complete parametrization of the real slice for generic torii only in the simplest case 
of a single spin. Fortunately, these real slices can be explicitely constructed for arbitrary n on 
level sets of the moment map containing very small strata, of minimal dimension m = (unstable 
critical points) or m = 1. The former case corresponds to isolated critical values and the later case 
to curves in the target space of the moment map. For these two situations, we have shown that it 
is possible to identify the physical strata of arbitrary dimensions in the corresponding level sets. In 
the case m = 0, the spectral polynomial has n + 1 double roots, and we recover the normal solitons 
constructed by Yuzbashyan [12]. When m = 1, Q2n+2(A) has only n double roots, but it defines 
a rational spectral curve, which allowed us to obtain the general form of the anomalous solitons 
discussed by Yuzbashyan [12] . We emphasize that these solitons can involve an arbitrary number 
of degrees of freedom, equal to the number of independent focus-focus singularites present in the 
quadratic normal form in the vicinity of the stratum of minimal dimension. 

This work raises several open questions. On the theoretical side, the most difficult seems to be 
to construct the real slice for level sets which contain a stratum of minimal dimension m > 2. For 
m = 2, the spectral curve can be uniformized by a non-singular curve of genus one, so we may hope, 
using Weierstrass functions, to find an explicit solution generalizing those we found for m = and 
m = 1. 

Another question is the quantum counterpart of these multi-mode solitons. If the quantum 
system is prepared in the separable state where each spin j is in the eigenstate of s? with the 
eigenvalue sej, ej = ±1, and the oscillator in the ground-state of b + b, how will such state evolve 
quantum-mechanically ? The problem is to decompose this initial state in the eigenvector basis 
of the quantum Hamiltonian. Some previous works have addressed this question in the case of 
a classical pinched torus of dimension 2 [17^ [T8l I19j . where the problem can be mapped into the 
evolution of a quantum wave packet for a single degree of freedom, prepared initially at the top of 
an unstable potential barrier. A complicated aperiodic sequence of pulses is obtained, each of them 
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being close to the classical monomode pulse of the n = 1 model. We suspect that the qualitative 
difference between classical and quantum evolutions will be more pronounced is the case of a critical 
level set of dimension larger than 2. Intuitively, we expect that quantum fluctuations perform some 
averaging over the configuration space of the multimode pulses whose result may be quite different 
from any classical pulse chosen in this ensemble. The theoretical challenge here is to extend the 
previous semi-classical analysis [TTJ, [181 CEH] to a system with several degrees of freedom. 

Finally, we hope that some of the aspects of the integrable dynamics discussed here will be 
evidenced experimentally. As already mentioned in the Introduction, one possibility is to consider 
cold gases of fermionic atoms in which an attractive interaction is switched on suddenly by sweeping 
the external magnetic field through a Feshbach resonance [HI [9] . The problem with this system is 
that it is most conveniently prepared in the ground-state of a weakly interacting Fermi liquid. 
After crossing the Feshbach resonance, this state becomes an unstable equilibrium point with only 
one focus-focus singularity in its normal form. So this gives rise only to the elementary single 
mode soliton. As shown in section El we need to start with an excited eigenstate of the diagonal 
part of the Hamiltonian to have a chance to observe multi-mode solitons. In this respect, another 
class of systems looks promising. Recently, several experimental groups have succeeded to couple 
a single quantum oscillator to a collection of electronic spins [20[ [21~1 [22"1 [23] . These experiments 
are in a regime where the spin ensemble behaves nearly as a macroscopic oscillator, so most-likely 
in a regime close to the ground-state of the coupled system. But it seems that such systems allow 
various manipulations |20[I25]. potentially useful for the long time storage of quantum information. 
We hope that they will provide a route to access also some of the interesting physics controlled by 
unstable critical points, which has been the main focus of the present work. 

9 Appendix: Normal form around critical points 

Critical points are equilibrium points for all the Hamiltonians Hj, j = 1, ■ ■ ■ , n + 1. At such points 
the derivatives with respect of all coordinates on phase space vanish. We have 2 n critical points 
located at: 

6 = 6 = 0, sf = 0, s] = ejS, e,- = ±1 (108) 

When we expand around a configuration given by eq. fl!08j) . all the quantities (6, b, sj~, s~) are 
first order, but s 2 - is second order because Sj = sej — ^s^s~ + • • • , ej = ±1. It is then simple 
to see that all first order terms in the expansions of the Hamiltonians Hj vanish. In particular, 
we can replace the dynamical generating function A(X) defined by eq. (J10[) by its non-dynamical 
approximation: 

a(A) = 2A + E _£^ 

j J 
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We can expand the Hamiltonians Hj around the equilibrium points eq. (|108l) and write them in 
normal form. Normal modes are obtained as follows. Consider the equation: 

a(E) = 2E + ^2 J^ = °i "Classical Bethe Equation" (109) 

3 J 

This is an equation of degree n + 1 for E. Calling Ei its solutions, we construct in this way n + 1 
variables C{ = C{Ei) and conjugated variables Bi = B(Ei). We have: 

{d,Cj} = 0, {B i ,B j } = (110) 

and 

{B i ,C j } = 2ia' i 8 tJ (111) 

where d i stands for a'(Ei). Up to normalisation, these are canonical coordinates. It is simple to 
express the quadratic Hamiltonians in theses coordinates: 

This has the correct analytical properties in A and in particular, there is no pole at A = Ej because 
a{Ej) = 0. Expanding around A = oo we get: 



H n+ i = s 22 e k + 2_J l^^iCi 

k 

and computing the residue at A = ej, we find: 



H 



Hj = sej 



^ + E 



se k 



r ^ - ek 



+ 7 — r —BjCj 

^2a' i e j -E i 



The physical Hamiltonian is then: 



H = H cp + J2 l ^^B i C i (113) 

i ' 

where H cp is the total energy at the critical point. This expression, together with Poisson brackets 
(IllOP and (jllip , shows that that eigenfrequencies for the linearized equations of motion under the 
Hamiltonian H are ±(w + 2Ei). 

We can invert these formulae: devide eq. (|112p by A — Ej and take the residue at A = Ej. Since 
a(Ej) = we get: 

hvL 2 (E J ) = B J C, 
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or explicitly: 

n 2H n s 2 

Bi c i = i Ei +i H n+1+ Y, w ^- + Y, JW — r2 

k=l ] K fc=l V 3 K ' 

Let us express now the first order quantities b, b, si", s~ in terms of the normal coordinates B 
C{. For this, we can reconstruct the generating functions B(X) and C(A). In fact: 

where N(X) is a polynomial of degree n and we know its values at the n + 1 classical Bethe roots: 

N{E i ) = B i J[(E i -e k ) 

k 

Hence we can reconstruct it by Lagrange interpolation formula: 






Now: 



so we can write: 



and similarly: 



, 2 !W*-%: 



n fe (^-^) 

rm _,Vr 1 IIfc^(A - E k ) 



a,- *-r a,- 



Computing the residue at A = ej we get: 






s - = 2 J^ ^ > ^ 



HMf( e J - efc ) V a i( e i - ^) 



We can simplify this formula. Using the definition of classical Bethe roots 



v A - e ^ n fc (A-e fc ) 



(114) 



(115) 



Computing the leading terms at oo we find: 

&=£-*«, 5=E^ ci (116) 



5G 



and calculating the residue at A = e,- we find: 



sen = 2 



3 ~ U&jfa ~ e *) 

hence: 

The reality conditions are easily expressed in this approximation. From B(X) = C(X*)* , we get: 

Bf = Cf 

where Cj stands for C(Ei). From Eq. f)115j) . it is possible to describe the pattern of separated 
variables \ a corresponding to the vicinity of the critical point on the critical torus. Using Eq. (|113|) . 
we see that this critical torus is defined, within the quadratic approximation, by: 

Bid = 

for any i. For a real root Ei, the reality condition implies Bi = C* so B{ = d = and we have one 
separated variable \ a frozen at Ei . Let us now consider a complex root Ei . Using the decomposition 
Eq. (jll3p of the physical Hamiltonian in normal modes, and Eqs. (|110|llip for Poisson brackets, we 
find: 

dtCi = {H, d} = i(uj + 2Ei)Ci 

so that: 

d(t) = Ci(0)e i( -" +2Ei)t 

The tangent cone to the critical torus at the critical point is composed of two hyperplanes. The first 
one corresponds to unstable directions, for which Ci(t) grows as time increases, that is QEi < 0. 
The second one is associated to stable directions, given by QEi > 0. Generically, on the unstable 
hyperplane, the expression ()115[) for C(A) is dominated at large times by the term such that QEi < 
with largest absolute value. Eq. (I115P shows that in this case, all the separated variables are located 
on the n roots Ej different from Ei. On a codimension one subspace on the unstable hyperplane, the 
coefficient Cj(0) of this leading eigenvalue vanishes, and the subleading root becomes non- frozen. 
More generally, C(A) is dominated at large times by the term for which Im Ei is negative and of 
maximal absolute value among the i's such that Cj(0) 7^ 0. 

Similarly, the stable hyperplane exhibits a nested sequence of linear subspaces, each of which 
corresponds to freezing the separated variables on all the roots Ej excepted one root Ei with positive 
imaginary part. 



57 



References 



[1] R. H. Dicke, Coherence in spontaneous radiation processes, Phys. Rev. 93, 99 (1954). 

[21 J. M. Raimond, M. Brune, and S. Haroche, Manipulating quantum entanglement with 
atoms and photons in a cavity, Rev. Mod. Phys. 73, 565-582, (2001). 

A. Blais, R.-S. Huang, A. Walraff, S. M. Girvin, ad R. J. Schoelkopf, Cavity quan- 
tum electrodynamics for superconducting electrical circuits: An architecture for quantum 
computation, Phys. Rev. A 69, 062320, (2004). 

[4] M. Gaudin, La Fonction d' Onde de Bethe. Masson, (1983). 

[5] E. Yuzbashyan, V. Kuznetsov, B. Altshuler, Integrable dynamics of coupled Fermi-Bose 
condensates. Phys. Rev. B 72 (2005), p. 144524. 

[6] E. Jaynes, F. Cummings, Proc. IEEE vol. 51 (1963) p. 89. 

[71 M. Brune, F. Schmidt-Kaler, A. Maali, J. Dreyer, E. Hagley, J. M. Raimond, and S. 
Haroche, Quantum Rabi Oscillation: A Direct Test of Field Quantization in a Cavity, 
Phys. Rev. Lett. 76, 1800-1803, (1996). 

R. A. Barankov, and L. S. Levitov, Atom-Molecule Coexistence and Collective Dynamics 
Near a Feshbach Resonance of Cold Fermions, Phys. Rev. Lett. 93, 130403, (2004) 

R. A. Barankov, L. S. Levitov, and B. Z. Spivak, Collective Rabi Oscillations and Solitons 
in a Time- Dependent BCS Pairing Problem, Phys. Rev. Lett. 93, 160401, (2004) 

[10] L. H. Eliasson Normal forms for Hamiltonian systems with Poisson commuting integrals 
- elliptic case Comment. Math. Helvetici 65 (1990) pp. 4-35. 

1 O. Babelon, B. Doucot, Classical Bethe Ansatz and Normal Forms in the Jaynes- 
Cummings Model. arXiv:1106.3274 

[12] E. A. Yuzbashyan, Normal and anomalous solitons in the theory of dynamical Cooper 
pairing, Phys. Rev. B 78, 184507, (2008). 

[13] B.A. Dubrovin, I.M. Krichever, S.P. Novikov, Integrable Systems I. Encyclopedia of 
Mathematical Sciences, Dynamical systems IV. Springer (1990) p. 173-281. 

[14] O. Babelon, D. Bernard, M. Talon, Introduction to Classical Integrable systems. 

Cambridge University Press (2003). 

[15] E. Sklyanin, Separation of variables in the Gaudin model. J. Soviet Math., Vol. 47, 
(1979) pp. 2473-2488. 



58 



[16] M. Audin, Spinning tops, Cambridge University Press 1996. 

[17] R. Bonifacio and R. Preparata, Coherent spontaneous emission, Phys. Rev. A 2, 336 
(1970). 

[18] O. Babelon, B. Doucot, L. Cantini, A semiclassical study of the Jaynes-Cummings 
model, J. Stat. Mech. (2009) P07011. 

[19] J. Keeling, Quantum corrections to the semiclassical collective dynamics in the Tavis- 
Cummings model, Phys. Rev. A 79, 053825, (2009). 

[20] Y. Kubo et al., Strong Coupling of a Spin Ensemble to a Superconducting Resonator, 
Phys. Rev. Lett. 105, 140502 (2010). 

[21] D. I. Schuster et al., High- Cooper ativity Coupling of Electron- Spin Ensembles to Super- 
conducting Cavities, Phys. Rev. Lett. 105, 140501 (2010). 

[22] R. Amsuss et al., Cavity QED with Magnetically Coupled Collective Spin States, Phys. 
Rev. Lett. 107, 060502, (2011). 

[23] P. Bushev et al., Ultralow-power spectroscopy of a rare-earth spin ensemble using a 
superconducting resonator, Phys. Rev. B 84, 060501(R) (2011) 

[24] Y. Kubo et al., Hybrid Quantum Circuit with a Superconducting Qubit Coupled to a 
Spin Ensemble, Phys. Rev. Lett. 107, 220501, (2011). 

[25] Hua Wu et al. , Storage of Multiple Coherent Microwave Excitations in an Electron Spin 
Ensemble, Phys. Rev. Lett. 105, 140503, (2010). 

[26] John Williamson, On the Algebraic Problem Concerning the Normal Forms of Linear 
Dynamical Systems. American Journal of Mathematics, Vol. 58 No. 1 (1936), pp. 141- 
163. 

[27] V. I. Arnold, Mathematical Methods of Classical Mechanics, Springer, New- York, 1997, 
Appendix 6. 

[28] M. F. Atiyah, Convexity and commuting Hamiltonians. Bull. London Math. Soc. 14(1) 
(1982) pp. 1-15. 

[29] V. Guillemin, S. Sternberg, Convexity properties of the momentum mapping. Invent. 
Math. 67(3) (1982) pp. 491-513. 

[30] San Vu Ngoc Moment polytopes for symplectic manifolds with monodromy. Adv. Math. 
208 (2007), pp. 909-934 



59 



[31] H. Duistermaat, On Global Action-Angle variables, Comm. Pure Appl. Math. 33 (1980) 
pp. 687-706. 

[32] Alvaro Pelayo, San Vu Ngoc. Hamiltonian dynamics and spectral theory for spin- 
oscillators. ArXiv 1005.0439. 

[33] I. Krichever, "Hessians" of integrals of the Korteweg-De Vries Equation and Perturba- 
tions of Finite- Zone Solutions. Soviet Math. Dokl. Vol. 27 (1983), No. 3, pp. 757-761. 

[34] M. Audin, Hamiltonian Monodromy via Picard-Lefschetz theory. Comm. Math. Phys. 
229 (2002) pp. 459-489. 

[35] M. Zou, Monodromy in two degrees of freedom integrable systems, J. Geom. Phys. 10, 
(1992) p. 37. 

[36] N. T. Zung, A note on focus-focus singularities, Lett. Math. Phys. 60, (2002), pp. 87-99. 

[37] N. T. Zung, Another note on focus-focus singularities, Diff. Geom. Appl. 7, (1997), p. 
123. 

[38] R. Cushman and J. J. Duistermat, N on- Hamiltonian monodromy, J. Diff. Eqs. 172, 
(2001) p. 42. 

[39] O. Babelon, M. Talon, Riemann surfaces, separation of variables and classical and quan- 
tum integrability. Phys. Lett. A. 312 (2003), pp. 71-77. 

[40] P. Griffiths and J. Harris, Principles of algebraic geometry, Wiley, New- York (1978), 
chapter 2. 

[41] O. Babelon, D. Talalaev, On the Bethe Ansatz for the Jaynes-Cummings-Gaudin model. 
J. Stat. Mech. (2007) P06013. 



60 



